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Chapter 1 



Introduction 



In ordinary stars, energy is lost through electromagnetic radiation from the surface. 
In order to accomplish this, a gradient in the temperature will be set up. The centre of 
the star will be relatively hot and the surface layers relatively cool. As energy is radiated 
away from the star, it shrinks just sufficient to provide the energy loss, and an equal 
amount gets added to the internally stored energy. 

The contraction of the star leads to a continual heating of the stellar interior. A time 
will come when the temperature at the centre will rise to ~ 10^ K. At this stage, a 
new source of energy will appear in the interior : fusion of hydrogen nuclei into helium 
nuclei. Nuclear fusion releases enormous amounts of energy, which can prevent the star's 
gravitational collapse. The star attains hydrostatic equilibrium, and remains stable at 
a point on main sequence in the Herztsprung-Russel diagram, determined largely by its 
mass. The duration of this phase is very long (biUions of years). That is why most stars 
that we see in the sky are main sequence stars. 

In the above process, eventually a time will come when the supply of hydrogen in the 
star gets depleted. This leads to the following sequence of events : (1) a drop in the energy 
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production, (2) collapse of the star and a consequent rise of temperature in the core of the 
star (which is now mostly helium) and (3) fusion of helium to carbon and oxygen nuclei. 
This will lead to a second relatively stable phase of the star's life. However, unlike the 
first stage, a qualitatively new phenomenon now takes place. The energy release during 
the second phase of core collapse will not only lead to core heating but, in addition, cause 
an expansion of the outer layers of the star. As a result, the star gets a bloated shape 
(expected radius ~ 100 million km), and the temperature of the outer layers drops. This 
is called a red giant. The helium burning phase is expected from theory to last several 
million years. 

The depletion of helium again leads to a core contraction and a consequent rise in 
temperature. The further evolution of the star depends on the initial mass of the star. 

According to stellar evolution theory, for stars whose initial mass is less than about 
6 Mq {Mq — solar mass), the rise in temperature is insufficient to start the next fusion 
cycle (carbon — > neon). So, the collapse continues, and the energy released eventually 
blows off the remaining outer layers of the star. This produces a remnant stellar core 
and an expanding envelope. Such objects are called planetary nebulae, and some of these 
have been observed in the sky. The duration of this evolutionary phase is relatively short, 
about tens of thousands of years. 

The remnant stellar core does not undergo gravitational contraction indefinitely. To 
understand this, one has to recall the Pauli exclusion principle in quantum mechanics. 
This principle states that no two fermions {e.g. electrons, protons, etc.) can possess 
identical quantum numbers (like spin, charge, angular momentum, etc.) at the same 
space-time point. 
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The remnant core contraction implies a high density situation for the materiaL That is, 
the constituents of the core are squeezed into smaller and smaller volumes of space. When 
the pressure brought on by the high density exceeds the electrostatic binding energy of 
the electrons to the atomic nuclei, the bound electrons will get detached from the nuclei 
and form a gas of free electrons. Because the electrons must obey the Pauli exclusion 
principle, the net result of the density squeeze is that the dense electron gas will behave 
as if it is a stiff ball of steel, i.e. it cannot be further compressed after a certain point. 
In other words, a dense electron gas will exert a pressure that will resist gravitational 
collapse. This pressure is not thermal but entirely quantum mechanical in origin, and is 
called degeneracy pressure. Configurations of the remnant core whose degeneracy pressure 
balances the gravitational attraction are called white dwarfs. At this configuration, the 
diameter is very small (to correspond to the requisite high density), about 1% of the solar 
diameter, i.e. comparable to the earth's diameter, and is very hot. From this point in 
time onwards, the white dwarfs gradually cool down. Because of the small surface area 
and the peculiarity of its composition, the time scale of cooling of a white dwarf is very 
large, several biUions of years. Typical densities of white dwarfs are (10^ — 10^) g cm'^. 

For stars with initial mass in excess of GM©, the core is sufficiently massive so that when 
it undergoes gravitational contraction at the end of helium burning stage, temperatures 
become so high that new fusion processes can occur. These produce heavier atomic nuclei 
: carbon, oxygen, neon, etc. ... , upto iron. By the time iron nuclei are produced, there is 
a huge build-up of Coulomb repulsive forces because of the presence of a large number of 
positively charged protons. Furthermore, the binding energy per particle is very high for 
iron nuclei, so that these nuclei are very stable. As a result, no further fusion reactions 
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are possible, and the core starts to collapse again. 

This time around, the gravitational collapse becomes too strong (because of higher 
mass), and will overcome the electron degeneracy pressure. Two final stages of this 
collapse are possible : neutron stars and black holes. 

According to detailed calculations, when the density of matter in the remnant core 
reaches about 10^^ g cm~^, the composition is substantial amount of free neutrons (that 
is, unbound to nuclei), in addition to nuclei and the dense ambient electron gas. These 
neutrons will not undergo beta decay because of phase space barrier brought on by the 
presence of dense surrounding electrons, whose Fermi momentum gets to be pretty high 
(> 20 MeV). Beyond a density of about 10^'' g cm~^, the nuclei 'dissolve' because 
of depletion of their proton concentration due to inverse beta decay (which becomes 
energetically favourable at high densities) , and the composition of the core becomes mostly 
neutrons, with a small admixture of protons and electrons. If the core is not too massive 
(calculations suggest an upper limit of 2-3 solar masses), the degeneracy pressure and 
the repulsive force of neutron matter can balance the gravitational attraction. A stable 
configuration is then possible, and it is called a neutron star. Because of the high densities 
involved, a neutron star is expected to be very compact by stellar standards, with radius 
about (10-15) km. Since degeneracy pressure plays an important role for the stability of 
white dwarfs and neutron stars, these are sometimes referred to as degenerate stars (in 
which all fusion reactions have stopped, unlike ordinary stars). 

If the mass of the collapsing core is more than 3Mq, the collapse will go unchecked. 
There is no physical mechanism known that can provide enough repulsive forces to halt 
the collapse. This will lead to a singularity situation, where the density will be extremely 
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high with a large surface gravity such that not even hght radiation can emerge from its 
surface. This is called a black hole. 

Immediately after the discovery of neutron in 1932, it was visualized by Landau that 
very compact stable objects, made up entirely of neutrons could be formed in stellar 
collapse. Baade and Zwicky |l|] in 1934 suggested that such dense objects could be the 
remnants formed in the aftermath of supernova explosions. The first quantitative theo- 
retical estimates for the bulk properties of neutron stars, such as mass and radius, were 
given by Oppenheimer and Volkoff 0. Although the existence of neutron stars were 
predicted in 1932 and the theoretical estimates of their mass and radius were done in 
1939 (Oppenheimer and Volkoff 0), it was not until the discovery of pulsars 0] in 1968, 
that the real astrophysical importance of the neutron star idea was established. It is now 
generally accepted that pulsars are rotating, magnetized neutron stars. In order to have a 
proper understanding of the structure and the dynamic of pulsars, it is necessary to study 
the structure and composition of neutron stars in detail. Besides, the maximum mass 
possible for stable neutron stars is important in another astrophysical situation, namely 
the identification of a possible black hole in compact binary systems. 

Typical neutron star mass and radius are expected to be about IMq {Mq = solar 
mass = 2 X 10^^ gm) and 10 km respectively. This means that the average density 
inside a neutron star is of the order po or higher, where po = nuclear matter density 
^ 2.8 X 10^^ g cm~^. Most of the matter inside a neutron star is expected to be at 
densities much above po- The behaviour of such matter is not well understood. From 
heavy-ion collision data at intermediate energies, one now hopes to derive reliable nuclear 
equation of state for nuclear matter for p = (2 — 3)po [i]- The interactions that will be 
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important beyond this density are not known, although several theoretical models have 
been proposed during the last twenty years. An additional difficulty is the many-body 
aspect of the problem. A proper way to incorporate the many-body correlation effect in 
high density strongly interacting matter is necessary, but there exists as yet no consensus 
on this subject. 

There is general acceptance that neutron star interiors can be divided into four distinct 
density regimes characterized by different compositions |^, ^] : 

1. The density at the surface made up basically of ^^Fe nuclei bound in a lattice im- 
mersed in a sea of relativistic nondegenerate electrons, which ranges from 7.86 g cm~^ 
to a density of ~ 10'' g cm~^. 

2. The second density region is between 10'' g cm~^ to 4.3 x 10^^ g cm~^. The latter 
corresponds to the neutron drip point, at which point onwards continuum neutron 
states start getting populated. 

3. The third density region starts at 4.3 x 10^^ g cm~^ and it continues all the way 
upto nuclear matter density of po = 2.8 x 10^^ g cm^^. In this density regime, 
the matter consist of neutron-rich nuclei, neutrons (expected to be superfluid) and 
small number of protons and electrons. 

The density regions (1), (2) and (3) together define the crust of the neutron star 
because the nuclei are expected to be in lattice structure. The crustal depth from 
the surface depends on the equation of state and roughly corresponds to about 
(10 - 15)% of radius. 
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4. The fourth density regime corresponds to p > po- The central densities can be 
^ 10 Po- 

Oppenheimer and Volkoff considered neutron star matter to be non-interacting 
neutron matter and obtained a value 0.7Mq as the maximum gravitational mass and a 
corresponding radius ~ 17 km. Today it is believed that neutron star interiors are made 
up asymmetric nuclear matter, with perhaps small admixture of pions and hyperons. Since 
neutron star matter contains matter at density higher than the nuclear matter density, one 
should take into account, as far as possible, the nuclear forces at short ranges. In the last 
twenty years, many attempts have been made to calculate the equation of state p = p{p), 
p = pressure of high density matter, taking into account the nuclear forces at short 
distances. Most such realistic calculations give the neutron star maximum gravitational 
mass to be between (1.8 — 2)Mq. 

Clearly, it is important that interactions among neutrons be included in any realistic 
calculation of neutron star structure. The equation of state where the density > 10 po 
is poorly understood due to uncertainty in nuclear interactions, and has been a focus of 
much theoretical research in recent years. In this thesis, we make a study of properties 
of dense nuclear matter using a field theoretical approach, which has gained increasing 
importance in the last few years. 

This thesis contains two parts. Part one (chapters 2-5) deals with the nuclear equa- 
tion of state based on relativistic mean field theory. In chapter 2, we briefly review the 
equations of state in both non-relativistic and relativistic approaches. In chapter 3, we 
make a detailed study of the chiral sigma model and use it to derive an equation of state 
of asymmetric nuclear matter. Recently, Datta and Alpar 0] showed that the Vela pulsar 
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glitch data suggest stiffer rather than softer equation of state for neutron star matter. 
The chiral sigma model seems to possess such a desirable property at high densities. We 
have extended the chiral sigma model calculation to finite temperatures in chapter 4. We 
also discuss the derivative scalar coupling model proposed by Zimanyi and Moszkowski 
1^ some years ago. In chapter 5, we examine in the context of our model (as discussed 
in chapter 3), the question of a possible phase transition from hadronic matter to quark 
matter at high densities. We discuss the formation of strangelets of large mass in quark 
matter at high densities based on non-relativistic treatment in chapter 6. The detection 
of strangelets may be the most unambiguous way to confirm the formation of quark-gluon 
plasma in heavy ion collision experiment. The study of this chapter gives basic idea on 
the strangeness contents in the high density matter such as core of the neutron stars. 

In the second part of the thesis, we consider certain astrophysical applications. These 
are the modelling of stable neutron star structure, radial oscillations of quark stars and 
cooling rates of neutron stars if these objects possess (m, d, s) quark matter in their cores. 
The structure of neutron stars, quark stars and quark star oscillations are discussed in 
chapter 7. In the chapter 8, we discuss the limitation of previously suggested neutrino 
emissivity formula for the quark matter in the neutron stars. In the last chapter, we 
summarize the highlights and the main points of this thesis. 
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Chapter 2 

The equation of state of matter in 
neutron star interior : a review of 
previous work 

2.1 Introduction 

The equation of state p{p) (p=pressure, p=total mass-energy density) above the 
nuclear matter density (p > 3 x 10^^ cm~^) plays a crucial role to determine the 
equilibrium structure of neutron stars. A considerable amount of work has already been 
done in the last two decades on this subject. In order to set a proper perspective for 
the discussion in the following chapters, we give here a brief review of previous work on 
equations of state of matter in neutron star interior. Upto about nuclear density the 
equation of state is reasonably well known, but the central density of neutron star can 
be almost an order of magnitude higher. In this regime, the physics is unclear. Below 
nuclear density the nuclear gas is dilute and one can use the perturbative methods, while at 
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much higher densities we would be in the asymptotically free region of quantum-chromo- 
dynamics. Neutron star matter is usually divided into four general density regions ||: 

1. Near the surface region, where density is upto 10^ g cm~^, a lattice of bare nuclei 
(mainly ^^Fe) is immersed in a gas of relativistic and degenerate electron gas. The 
equation of state in this region is influenced by temperature and magnetic fields . 
But it hardly matters to the structure of neutron stars. 

2. The next region is the neutron-rich nuclei upto ^^^Kr, where the protons inside the 
nuclei undergo inverse beta decay. Due to the inverse beta decay, the electron gas 
occupies the lowest Fermi energy level, which allows the system to be in a lower 
value for the ground state energy. 

3. The third density region begins at about 4.3 x 10^^ g cm~^ is called the neutron 
drip point. In this region, some of the neutrons in the nuclei get detached from the 
parent nucleus. These neutrons are unbound and stable. 

4. In the final region, (density > 2.8 x 10^^ g cm~^), the individual nuclei merge into 
each other. So, the composition of the matter is expected to be a fluid of almost 
uniform neutron matter together with protons, electrons and possibly muons, pions, 
hyperons etc. Above this density, there is possibility of a phase transition from 
neutron matter to {u, d, s) quark matter. 

Because of large neutrino emissivity, neutron star matter is cold and degenerate. 
Therefore, we do not have to worry about temperature in deriving equation of state. 

The first three of the above density regions are reasonably well-understood. Compre- 
hensive account of the physics of these regions are given by Canuto |^ and Baym and 
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Pethick 0. All versions of equation of state of the fourth region (which is above the nu- 
clear matter densities) differ from each other, due to the inadequate knowledge of nuclear 
interactions as well as the lack of proper many-body techniques, that will be relevant to 
describe neutron star matter. In what follows, we briefly review the two major techniques 
that have been used : 

1. non-relativistic theories, 

2. relativistic field theoretical model. 



2.2 Non-relativistic approach 

In this approach one uses two-body potentials which are fitted to the nucleon-nucleon 
scattering data, as well as a three-body term. The form of the potential is chosen whose 
parameters are determined using the data of few-body nuclei and saturation properties 
of nuclear matter. For a non-relativistic model, the starting point is the two-body poten- 
tial. For the many-body method, there are two approaches, Brueckner-Bethe-Goldstone 
theory and the variational method |jTO|. Since the major constituents of matter inside 



the neutron stars are neutrons, we concentrate our attention on neutrons and their in- 
teractions. The experimental information regarding nucleon-nucleon scattering data and 
known properties of deuteron do not uniquely determine the nucleon-nucleon potential. 
Hence, for setting up the equation of state, it is required to have fit with known properties 
of equilibrium nuclear matter at saturation density such as binding energy per nucleon, 
compression modulus etc. Before we discuss the recent non-relativistic equation of state 
given by Wiringa et al. |Tl| that provides a reasonably good description, we briefly review 
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the historical equations of state at high densities in the following sections. 

/. Reid Model: 

The Reid model |T2[ is based on phenomenological nucleon-nucleon potentials, which 



fits the scattering data very well. It has been used extensively in the calculation of 
neutron star structure (reviews by Baym and Pethick ||l3l; Canuto fl^)- The calculation 



of equation of state p{p) using the lowest order constrained variational method seems to 
be accurate enough for the central part of the Reid potential (Pandharipande and Bethe 
T5|). But Pandharipande and Wiringa |T6[ calculated the nuclear matter properties using 



Reid potential and found that both the equilibrium density as well as binding energy were 
too large. Hence, this model is now considered to be unrealistic. 

//. Bethe- Johnson model: 

The short-range interaction between nucleons is not uniquely determined by the nucleon- 
nucleon scattering data. So, Bethe and Johnson [|1^ proposed a phenomenological po- 



tential model where they suggested several different potential models for nucleon-nucleon 
interaction, assuming various plausible strengths for short-range repulsion by fitting the 
scattering data. The maximum mass calculations for neutron stars with this equation 



of state are given in Malone et al. |T^. At high densities the constituents of matter are 
made up of nucleons (iV) and hyperons (Y). However, the hyperonic interactions are not 
included in full details. 

These two models (I and II) are phenomenological density independent static potential 
models. According to these, the nuclear matter energy at high densities increases linearly 
with the density. 
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///. Tensor Interaction model: 

The attraction between nucleons which comes from higher order contribution of the 
pion exchange tensor interaction was studied by Pandharipande and Smith [0. Their 
work is a generahzation of various tensor interaction models proposed earher by Green and 
Haapakoski ||20| . These interactions fit only the s-wave scattering data, and differ mainly 



in the strength of the short-range repulsion for which a specific form is presumed. The 
work by Smith and Pandharipande in pT| suggests that the low energy nucleon-nucleon 



scattering data can be explained by attributing all the attraction between nucleons to 
tensor interaction. It has been seen that the calculation with tensor interaction using 



lowest order variational and Brueckner methods (Green and Niskanen p2|) satisfies only 
half of the nuclear matter binding energy at saturation density. Thus the tensor interaction 
model can not explain the nuclear matter properties in a satisfactory way. 

IV. The mean field model: 



Pandharipande and Smith |^ assumed a model, called the mean field model, which 



states that the attraction between nucleons is due to the exchange of an effective scalar 
meson. The quadrupole moment of the deuteron as well as the phase shift in the ^Pq, ^Pi 
and ^P2 channels clearly indicate the presence of the one-pion exchange tensor force 
between the nucleons. However, a detailed analysis of the attractive interaction due to 
all possible tensor potentials (Smith and Pandharipande [^]) suggest that it is almost 
independent of the spin and isospin of the interacting nucleons, and thus its contribution 
in matter could be similar to that due to coupling of nucleons to a scalar field. The 



scalar field, treated in the mean field approximation, was used by Walecka |2|], which 
we shall discuss later in detail. Following the Walecka model, the nucleons moving in a 
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mean scalar field are assumed to interact by a central potential generated by uo, p and vr 
mesons exchanges. The potential approximation may be appropriate for u and p vector 
fields, while the central parts of the pion exchange potentials have a negligibly small effect. 
The mean field approximation is unjustified for the vector field because the u-p exchange 
potential have a range ~ 0.2 fm, which is much smaller than the mean interparticle 
spacing ~ 1.2 fm. The short-range corrections induced by the u-p exchange potentials 
are treated by variational method using a hypernetted chain formalism (Pandharipande 
and Bethe |2^). The coupling constant of u, p and a are calculated from the nuclear 
matter binding, symmetry energy and the equilibrium density. The incompressibility 
parameter obtained by this theory is ~ 310 MeV. The nucleon-nucleon scattering data 
can not be explained by an interaction which has been used in the mean field model. But 
this model satisfies all the empirically computed properties of nuclear matter. 

The attraction between nucleons in the above two models, i.e., tensor and mean field 
decrease with increasing density. However, this is the general characteristic of microscopic 
models, which are based on the mean field theoretic calculations. The draw back of these 
two models is that at small densities the energy is proportional to the density, whereas at 
higher densities it tends to saturate. 

V. Friedman and Pandharipande model: 

Friedman and Pandharipande PB[ gave a model for the equation of state of dense 



neutron and nuclear matter, where they used the variational method as suggested earlier 
by Pandharipande [110] to calculate the equation of state for a wide range of density. 



In their model, they used improved phenomenological nucleon-nucleon interactions 127 



28| . The two nucleon interaction has short-range and intermediate range parts, and also 
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possesses the pion-exchange contribution. There is also a contribution due to a three- 
nucleon interaction which has a comphcated form (it is a function of strength parameters, 
the interparticle distance and also the alignment angles). The parameters involved in the 
three-nucleon interactions are determined by reproducing the equilibrium density, energy 
and incompressibility of nuclear matter based on variational calculations. This model fits 
well the nucleon-nucleon scattering cross-section data, the deuteron properties and the 
nuclear matter properties. 

VL Wiringa, Fiks and Fabrocini model: 

Wiringa et al. [|ll| proposed a model which is most firmly based on available nuclear 



data. This model improves on the earlier work by Friedman and Pandharipande |^6[. In 
this model the two-nucleon potential is taken to be the Argonne vl4 (AV14) or Urbana 
(UV14) potential. Both have the identical structure, but differ in the strength of the 
short-range tensor force. They are called vl4 models because they are the sum of vl4 
operator components (like aj.cTj, TTj.TTj, etc.). Each of these components has three radial 
pieces which includes the long-range one-pion exchange, an intermediate-range part that 
comes from the two-pion exchange processes, and a short-range part, coming either from 
the exchange of heavier mesons or overlaps of composite quark systems. All the free 
parameters are fitted to nucleon-nucleon scattering data and deuteron properties. 

For the three-nucleon interaction, the Urbana VII potential is used, which has a two- 
pion exchange part and intermediate-range repulsive contribution. Calculations have been 
carried out for the Urbana f 14 plus three-nucleon interaction (TNI) model of Lagaris and 
Pandharipande ||29 . 



The many-body calculations are based on the variational principle where one uses 
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the technique called the Fermi hypernetted chain-single-operator chain (FHNC — SOC) 



integral equations ^ 



The authors obtained the nuclear matter saturation properties such as binding energy, 
saturation density and incompressibility parameter value of nuclear matter at saturation 
density are -16.6 MeV , 0.157 fm'^ and 261.0 MeV respectively for UV14 plus TNI 
model. This is a reasonable improvement over all the other non-relativistic approaches. 

In their study, one notices that the sound velocity s in the medium given by the equa- 
tion of state (in parametric form) for beta-stable (n, p, e i.e., no hyperons) matter based 
on non-relativistic approach violates causality above p = lfm~^. This is an undesirable 
feature of this method at high densities. However, they predicted the maximum neutron 
star mass to be 2.2Mq and for I.AMq neutron star, the central density turns out to be 
substantially below lfm~^ which is quite realistic. It may be noted that the modern po- 
tentials supplemented with reasonable three-body interactions yield very similar models of 
neutron star structure parameters. For the 1.4Mq models, the radius is ~ 10.4 — 11.2 km, 
and the central density is about 6po- 



2.3 The relativistic approach 

The shape of baryonic potential is not known at very small interparticle separations 
(< 0.5/m). Also, it is not clear that the potential description will continue to remain 
valid at such short ranges. Moreover, in the neutron star interiors, the Fermi momentum 
of the degenerate neutron is large. Therefore, the non-relativistic approach may not be 
adequate. In the recent times, the relativistic approach has drawn considerable attention. 
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In the relativistic approach, one usually starts from a local, renormalizable field theory 
with baryon and explicit meson degrees of freedom. These models have the advantage of 
being relativistic, however, one drawback is that one does not know how to relate it to 
nucleon-nucleon scattering data. The theory is chosen to be renormalizable in order to fix 
the coupling constants and the mass parameters by empirical properties of nuclear matter 
at saturation ( binding, density, compression modulus, effective mass and symmetry en- 
ergy). As a starting point, one chooses the mean field approximation (MFA) which should 
be reasonably good at very high densities (a few times nuclear) [^. In the second step, 
one includes one-loop vacuum fiuctuations which leads to what is called the relativistic 
Hartree approximation (RHA) . This approach is currently used as a reasonable way 
of parametrizing the equation of state. 

The main features of this approach can already be seen in a simple model. Assuming 
that a neutral scalar meson field (0) and a neutral vector field (V^) couple to the baryon 
current by interaction terms of the form 




the Lagrangian can be written as {h 



1 



C = i>[r{^^,~g,V,)-{M 



(2.2) 



where 



F, 



(2.3) 



From Euler-Lagrange equation 
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d , dC . dC , 
[^T^TT^] - ^ = (2-4) 



dxf^ d{dqi/ dxf^ ) dq, 
where Qi is one of the generahzed coordinates, one yields the field equations 

id^d^ + ms')<P = gsipip, 

[Yitdf, - gvV^) - (M - g,<P)]'iP = 0. (2.5) 

The MFA consists in replacing the meson-field operators by their expectation values : 

<0>=0„ <V^>=5^oVo, (2.6) 

if we consider a static uniform system. Then, the meson field equations become 

0, = ^ < >= ^ps, (2.7) 
rrig rrig 

K = 4 < >^ ^PB. (2.8) 

In this approximation, the nucleon field operator satisfies a linear equation : 

Hrdi, + g,YVo + M*]ij = 0, (2.9) 

where M* — M — gs(j)o is the effective mass of the nucleon and the Lagrangian density 
takes the form 

CM:Fr = i^ih^^d^ - gvYVo - M*]^ - ^my, + \myi (2.10) 

The quantization here in this case is straightforward. In the mean field approximation, 
the conserved energy-momentum tensor is 
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T,. - -g,XMrr + 9^^^^^^- (2-11) 

Assuming that the bulk neutron star matter can be treated as a uniform and perfect fluid, 
we can write the energy density of the system as 

e =< T,, > (2.12) 

and the pressure as 

p^^<Tii>. (2.13) 

The ground state of nuclear matter is characterised by Fermi momentum kp which is 
related to the baryon density ps s& : 



Here 7 = 4 for symmetric nucleon matter, 7 = 2 for neutron matter. The mean field 
00 (or M*) is determined self-consistently as 

M*^M-g,<f>,^M-^p,^M-^-L- d'k---, (2.15) 



m1{2TiY Jo E*{k) 

where E*{k) — (/c^ + M*^)^/^. This leads to a transcendental equation for a given kpipB)- 
It turns out that only the ratio g1/m1 and g1/m1 enter in the equation of state. These 
parameters are fixed such that the nuclear matter will give rise to the correct binding 
energy and saturation density : 

E - BM 

( )o = -15.75 MeV, (2.16) 

rS 

k"p = 1.42 fm-^ (7 = 4). (2.17) 
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Thus 



Cl^gl{—) = 2mX (2.18) 



Cl = gl{—) = 195.9. (2.19) 

ml 



With these values one can compute the equation of state for neutron matter. In this 
simple model there is a first order phase transition (similar to the liquid-gas transition in 
the van der Waals' equation of state). At asymptotically high densities, the velocity of 
sound approaches the velocity of light. 

One can make the model more realistic by including for instance, the interaction of p- 
mesons (charged vector mesons). In such model (QHD-II of Serot and Walecka [^) the 
p meson stiffens the equation of state at relatively low density and causes the gas-liquid 
phase transition to disappear. 

The above field theoretical treatment is now-a-days referred to as the Walecka model 
2^ . In 1974, Walecka first proposed this mean field model, where he chose the coupling 



constants in such a way that it fitted the nuclear matter binding energy and saturation 
density. In this model, the value of nuclear matter incompressibility at saturation is quite 
high. The isospin triplet vector meson p, was not included, but is of relevance in neutron 
star interior (ra, p, e) matter. The extension of this model will be discussed in the next 
section. 

/. Glendenning Model: 



Glendenning [31] presented a relativistic field theoretical model for densities near as 
well as above the nuclear matter density. This includes isospin-asymmetric baryon matter. 
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In this model was included the nucleons, the mesons : tu, p and vr, and /i^ particles 
along with the self-interaction of the scalar meson field a. The self-interaction of a-field 
was chosen to be of the form : 



U{a) = {Kri^ + -^g,a){g„af (2.20) 



where m„ is the nucleon rest mass and 6, c are parameters. The coupling constants and 
the parameters 6, c were determined by satisfying the following bulk properties of nuclear 
matter i.e., 

saturation particle density = 0.145 firr^ 

saturation binding energy = -15.96 MeV /nucleon 

asymmetry energy = 37 MeV 

and the incompressibility = 280 MeV 

The author chose two different sets of parameters, one of which gives a soft and other 
a stiff equation of state at high densities, for fits to nuclear physics. 



In 1985, Glendenning extended the above formalism to derive an equation of state 
of hyperonic matter. He used the same constraints as earlier to fix the various parameters. 
Here, he included the K and K* meson exchanges and also the self-interaction form as 
earlier. The coupling strengths of the isobar A and the hyperons (A, S^'", S^'*^) to the 
mesons are taken from the work of Moszkowski [Rll . This work indicates that cores of the 



heavier neutron stars are dominated by hyperons and the total hyperon population for 
such stars are 15% — 20%, depending on whether pions condense or not. The advantage of 
this model is that one can make the equation of state soft or stiff depending on the input 
parameters, namely the incompressibility value. In this case the following four separate 
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hyperonic equation of state matter are considered. 
Case 1 : n, p, hyperons, A, e", vr"; 
Case 2 : n, p, hyperons, A, e^, /i^; 

Case 3 : like case 2 but with universal coupling of hyperons; 
Case 4 : hke case 3 but no p- exchange {g^ = 0). 

We describe now the results of a more recent analysis for the following (cr, u, p) 
model, whose Lagrangian is given as 

B ^ 

1 11 

-\p,ur + I^Ip.-P" - \bmn{g.<yf - (2.21) 

Here the sum runs over the baryons A^, A, S, S, A, etc. 

For uniform matter, the theory depends only on the ratios ga/ma, g^^I'm^, 9p/^p and 
the self-interaction coefficients 6, c. These five quantities are determined empirically, in 
both the MFA and the RHA ( Table 1 in Ref. |^). From this calculation one notices that 
the equation of state for stable neutron star matter will be softer than the pure neutron 
matter. This softening is the result of the replacement of energetic neutrons by hyperons 
at rest. It is most remarkable that the maximum allowed neutron star mass depends very 
sensitively on this coupling. This point has been further investigated in Ref. with the 



conclusion that even if we have a perfect knowledge of the nuclear equation of state upto 
about 2po, there is still a large uncertainty in the maximum neutron star mass because 
of our lack of knowledge of the hyperon interactions. This is one of the most relevant 
"interface" of astrophysics and nuclear physics. However, further laboratory experiments 

22 



are needed to reduce the large uncertainties. 

//. Alonso and Cabanell model: 



Alonso and Cabanell fS^] gave a model in which the nucleons interact via scalar (cr) 
mesons, pions (vr), and vector {uj and p) mesons. The model is solved in the renormalized 
Hartree approximation. The authors give two sets of equation of state (I, II). 

The equation of state I is derived by fitting the properties of symmetric nuclear matter 
at nuclear density in a satisfactory way. This is in good agreement with the equation of 



state obtained by Baym, Bethe and Pethick [38| in the region above the neutron drip 
for neutron matter. The value of the nuclear incompressibility obtained in this model is 
too high (460 MeV). This large incompressibility is perhaps due to the absence of a- 
self-interaction term in the Lagrangian. 

The equation of state II is derived using the chirally invariant a- model Lagrangian, 
coupled to uj and p mesons and having an explicit symmetry breaking term. This equation 
of state fits accurately all the known properties of symmetry energy of nuclear matter at 
nuclear density. The value of the nuclear incompressibility at saturation density comes 
out to be 225 MeV. 

For these two equations of state, the total number of adjustable parameters are rather 
large (12 for set I and 15 for set II). A peculiar feature of these models is that both (I, 
II) show a dip in equation of state. These dips do not represent any phase transition for 
the system. 

///. Chiral sigma model (Glendenning): 

Chiral symmetry is a good hadron symmetry, ranking only below the isotopic spin 
symmetry |^9[|. For this reason, it is desirable to have chiral symmetry in any theory of 
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dense hadronic matter. At the same time, the theory should be capable of describing the 
bulk properties of nuclear matter. But a theory that satisfies both the conditions is not 
available. 



In 1986, Glendenning |^ derived the equation of state based on the mean field ap- 
proximation. Here, the equation of motion for the mean fields is also derived in the 
general case for the chiral sigma model supplemented by a gauge massless vector meson 
in interaction with the other hadrons, including baryon resonances. Here, omega meson 
is considered with dynamical masses. In this theory, the gauge field uj^ of massless vec- 
tor meson is introduced into the chiral sigma model through the covariant derivatives. 
Moreover, the linear term of a field in it is the symmetry breaking term, by which the 
pion acquires a finite mass. In addition to the a;-mesons, the theory has scalar meson 
and the pseudoscalar pion. This theory has two parameters, which are determined from 
the saturation density and binding energy per nucleon in normal symmetric matter. The 
value of the incompressibility is rather very large i.e., K = 650 MeV. 



In a subsequent paper, Glendenning [41| extended the chiral sigma model based on 



the mean field approximation where the a;-meson does not have a dynamically generated 
mass, even if one considers the vacuum renormalization correction to the theory. In this 
model, he considered only the normal non-pion condensed state of matter and included 
hyperons in beta equilibrium with nucleons and leptons. He fitted the parameters of the 
theory to obtain two equation of state corresponding to two compression moduli, a "stiff' 
{K = 300 MeV) and a "soft" {K = 200 MeV), by reproducing correct nuclear matter 
properties. 
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IV. Chiral sigma model (Prakash and Ainsworth): 

Prakash and Ainsworth proposed an equation of state based on the chiral sigma 
modeL They examined the role of the many-body effects provided by the chiral sigma 
model in the equation of state of symmetric nuclear matter and neutron rich matter. 
They include the a- meson one-loop contributions, but the isoscalar vector field is not 
generated dynamically, so that its role is reduced to an empirical one. A set of equation 
of state is constructed, each of which fits the empirical saturation density, the binding 
energy and the symmetry energy of nuclear matter. They obtained a value of compression 
modulus at saturation nuclear matter density for symmetric matter that is different from 
the experimental expected value. So, they allowed a variation in the values of coupling 
constants arbitrarily so that the theory will give rise to the desired value of compression 
modulus. The vector field plays no role in determining the value of the effective mass of 
the nucleon in such an approach. 

V. Baron, Cooperstein and Kahana model: 

Baron et al. gave a phenomenological model for neutron matter equation of state. 
The form of the nuclear partial pressure is 



P = ^^[u^-il (2.22) 



where the baryon density compression factor u = p/po, measured with respect to the em- 
pirical symmetric nuclear matter saturation density, po = 2.8 x 10^^ g cm~^ = 0.16 fm~^. 
Here, Kq is the nuclear incompressibility at saturation, i.e. u = 1, while 7 is the extremely 
high-density adiabatic index. 
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Later, Cooperstein modified the phenomenological equation of state to tlie follow- 
ing form 

KnPo 



97 



-M^, (2.23) 



P 

e = po{m„ + mi)u H -, (2.24) 

7-1 

where the integration constant mi is given by 

'^l — ~l" Egym 7^ 7T, (2.25) 

97(7 - 1) 

with the empirical value Enm = —IQMeV and for definiteness E^y^ = 36 MeV. Here 
the factor u'^ is replaced by m'*' — 1 for nuclear matter saturation. The neutron matter 
compression modulus is denoted by Kn and m„ is the nucleon mass. The value of 
is not interpreted as the correct saturation density but is represented as a high density 
neutron matter parameter. This is connected with Kq through the density dependence 
of the symmetry energy. One can find different equations of state by varying the values 
of K]\f and 7 in this theory. For particular values of and 7, these phenomenological 
equations of state fit well the non-relativistic equation of states, such as that of the Bethe 
and Johnson or Friedman and Pandharipande ||26| . 

VI. Rosenhauer et al. model: 



Rosenhauer et al. [45] refer two commonly used parameterizations for the hadronic 



equation of state, namely that of Sierk and Nix |^ 

2K 



esNin) = ^{^n/no-lf (2.26) 



and the quadratic form is 
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^qH = Y^(^-^o)W- (2.27) 

The symbol esN^n) and eq^n) represent energy density of the two equations of state, 
namely, Sierk and Nix as well as quadratic form respectively. Here, hq is the normal 
nuclear matter density and K, the compression modulus characterizing the properties of 
nuclear matter at densities n > uq. In these two parametric equations, ecomp{n) {i.e., 
£sN{n) or EQ^n)), the large K value indicates the stronger repulsive nature of nucleon- 
nucleon interaction. 

The total energy density is then expressed as 

e{n) = n[ecomp{n) + Wo + mn + W,yms], (2.28) 

where the binding energy per nucleon at normal nuclear matter density, (no = 0.145 fm~^) 
is Wo = —16 MeV, rrin being the rest mass of the neutron and Wgym = 32 MeV, 
the symmetry energy of neutron matter estimated from the liquid drop model. The 
symmetry energy Wsym and the binding energy Wq determine the properties of matter 
at saturation within this phenomenological approach. However, Ecompin) incorporates all 
density dependent effects. 

VII. Zimanyi and Moszkowski model: 

Zimanyi and Moszkowski proposed a model similar to the Walecka model but 
using the derivative scalar coupling (DSC) for the scalar field. The interesting feature of 
this model is that the equation of motion of the scalar field becomes non-linear without 
introducing any extra parameter. This model gives rise to a reasonable value of the 
effective nucleon mass and satisfactory value for the incompressibility of nuclear matter 
at saturation density. 
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Glendenning et al. extended the above model to hyperonic matter and applied to 
compute a number of neutron star properties. In place of the purely derivative coupling 
of the scalar field to the baryons and vector meson fields of the above (DSC) model, they 
here coupled it by both Yukawa point and derivative coupling to baryons and both vector 
fields. This improves the value of the compression modulus and effective nucleon mass 
at saturation density compared to that obtained in earlier calculation. Also, they include 
the fh meson to account for the asymmetry effect. 
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Chapter 3 



Chiral sigma model 



3.1 Introduction 

The chiral symmetry is a good hadron sjTiimetry, which ranks only below the isotopic 
spin symmetry Due to this reason, it is expected that any theory of dense matter 
should possess it. If one considers the chiral symmetry in dense matter theory, it should 
be capable of describing the bulk properties of nuclear matter such as binding energy per 
nucleon, saturation density, compression modulus and symmetry energy. So far, there 
has not been any theory possessing chiral symmetry and describing all the nuclear matter 
properties. In recent years, the importance of the three-body forces in the equation of state 
at high densities has been emphasized by several authors [^, This gives theoretical 
impetus to study the chiral sigma model, because the non-linear terms in the chiral sigma 
Lagrangian can give rise to the three-body forces. 

A chiral Lagrangian using the scalar field (the so-called sigma model) was originally 



introduced by Gell-Mann and Levy |51| as an example to illustrate chiral symmetry and 
partial conservation of axial current. The importance of chiral symmetry in the study of 
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nuclear matter properties was emphasized by Lee and Wick [Q. The non-hnear terms 
of the chiral Lagrangian can provide the three-body forces, important at high densities 
(Jackson, Rho and Krotscheck [^]; Ainsworth et al. [Q), and can be relevant in appli- 
cations to neutron star structure and supernova collapse dynamics. 

The usual theory of pions leads to a theory of nuclear matter that does not possess 
the empirically desirable saturation property. For this reason, the isoscalar vector field 
is introduced in the theory via the Higgs mechanism. This way it becomes possible to 
have a saturating nuclear matter equation of state (Boguta With the availability 

of experimental estimate for the incompressibility parameter of nuclear matter (denoted 
by K), there have been attempts to reproduce the desirable value of K (about 200 — 300 
MeV) using the sigma model. In the 'standard' sigma model, the value of K turns out to 
be quite large, several times the above-mentioned value for plausible values of the coupling 
constants involved, and can be reduced only by introducing in the theory, terms due to 
the scalar field self-interactions and/or vacuum fluctuations with adjustable coefficients. 

There have been several earlier papers that have employed the chiral sigma model to 



obtain the equation of state of high density matter. Glendenning derives the chiral 
sigma model equation of state with normal nuclear matter saturation and then applies to 
neutron star structure calculations. He then studies the finite temperature behaviour of 
this model. A liquid-gas phase transition occurs in this model below T ~ 23 MeV . The 
equation of motion is derived in the mean field approximation, supplemented by a gauge 
massless vector meson in interaction with the other hadrons, including baryon resonances. 
The gauge field of a massless vector mesons is introduced through the covariant derivatives 
into the chiral sigma model. The linear term of sigma field in the Lagrangian is the 
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symmetry breaking term by which the pion acquires a finite mass. There are scalar meson 
and pseudoscalar pion in the theory, in addition to uo- meson, which contribute a repulsive 
force and carry spin and isospin. Also, he considers this model with dynamical masses but 
not with the p- meson and its isospin symmetry infiuences. The required parameters in 
this theory are calculated from the saturation density and the binding energy per nucleon 
in normal symmetric matter. The value of the compression modulus obtained from this 
model is too large (650 MeV). 



In another paper, Glendenning |^ extends the chiral sigma model to normal nuclear 
matter and neutron star matter by considering the vacuum renormalisation corrections. 
The equation of state is derived by including p- mesons and hyperons in equilibrium with 
nucleons and electrons and muons. In this calculation the uj- meson does not have a dy- 
namically generated mass. He has considered two sets of coupling constants corresponding 
to two compression moduli, a 'stiff' {K = 300 MeV) and 'soft' {K = 200 MeV) equation 
of states for nuclear matter, that yield the empirical saturation density, the binding as 
well as the symmetry energy. In both cases, the hyperons influence on the equation of 
state substantially. 

An equation of state based on the chiral sigma model is also considered by Prakash 



and Ainsworth [42|. They examine the role of the many-body effects provided by the 



chiral sigma model in the equation of state of symmetric nuclear matter and neutron-rich 



matter. Previously, Matsui and Serot considered the role of only the nucleon vacuum 
fluctuation terms at the one loop level. But here, these authors consider both the nucleon 
as well as the a meson vacuum fluctuation terms at one loop level in the chiral sigma 
model. As a result, the meson loop shifts the saturation density and thus incompressibility 
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increases. In this theory, the isoscalar vector field is not generated dynamically. Its role is 
an empirical one. They have constructed a family of equation of states, each of which fits 
the empirical saturation density, the binding energy and the symmetry energy of nuclear 
matter. They allow variation of the coupling constants arbitrarily in their calculations 
by fitting the other parameters from nuclear saturation properties, so that it is possible 
to obtain any desired value of the nuclear matter compression modulus. Finally, they 
calculate neutron star structure using the neutron matter equation of state for p, n and 
e systems. 



We have developed the SU (2) x SU{2) chiral sigma model to describe nuclear 
matter and neutron star matter. We have adopted an approach in which the isoscalar 
vector field is generated dynamically. Inclusion of such a field is necessary to ensure the 
saturation property of nuclear matter. The effective mass of the nucleon thus acquires a 
density dependence on both the scalar as well as the vector fields, and must be obtained 
self-consistently. We do this using the mean-field theory where all the meson fields are 
replaced by their uniform, expectation values. To describe the nuclear matter we have two 
parameters in the theory : (i) the ratio of the coupling constant to the mass of the scalar 
and (ii) to the isoscalar vector fields. This procedure also gives a relatively high value for K 
at the saturation density. Although this is an undesirable feature as far as nuclear matter 
at saturation density is concerned, it need not be viewed as a crucial shortcoming for our 
purpose here in view of the fact that a fit to at saturation does not tell us what the slope 
of the equation of state should be at densities > {ug = saturation density) (Prakash 



and Ainsworth [|^; Horowitz and Serot [Q; Stock [^]; Baym Ellis, Kapusta and 



Olive 15^). For neutron star structure, in which we are most interested as an application 
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of our equation of state, this high density regime plays the most important part. The 
important feature of our work is that the lo- meson mass and couphng constant are not 
treated empirically, but in a fully self-consistent manner. To describe neutron star matter, 
we include the coupling to the isovector p- meson, the coupling strength being determined 
by requiring a fit to the empirical value of the symmetry energy. 

3.2 The Model 

The Lagrangian for an SU{2) x SU{2) chiral sigma model that includes (dynami- 
cally) an isoscalar vector field {uj/j) is (we choose h — 1 — c) : 

C = ^{d^^.d'^^ + d^ad^a) - ^(V.^ + (7^ - xlf 

11 2 

^Qaipia + z-zsT". + - g^'^^.uj'')^), (3.1) 

where F^^y = d^uj^, — d^uj^, ip is the nucleon isospin doublet, "t? is the pseudoscalar pion 
field and a is the scalar field. The vector field lO/j, couples to the conserved baryonic current 
Jn — '4^1n'4^- The expectation value < jo > is identifiable as the nucleon number density, 
which we denote by p^- 

The interactions of the scalar and the pseudoscalar mesons with the vector boson 
generates a mass for the latter spontaneously by the Higgs mechanism. The masses for 
the nucleon, the scalar meson and the vector meson are respectively given by 

M = g^Xo] 
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(3.2) 



where Xq is the vacuum expectation value of the sigma field, A = (rrio-^ — '^7r^)/(2/jr^) 
(m7r = pion mass and /tt is the pion decay coupling constant) and and are the 
coupling constants for the vector and scalar fields respectively. 

To derive the thermodynamic quantities of the system of degenerate nucleons, char- 
acterized by the nucleon number density {pb) or equivalently the Fermi momentum 
kp = {Qt^"^ Pb hy^^{l = nucleon spin degeneracy factor), we need to know the depen- 
dence of the meson fields on p^. For this, we resort to the mean-field approximation. 
This approach has been extensively used to obtain field theoretical equation of state mod- 
els for high density matter. In this approximation, expected to be valid for degenerate 
matter at high densities, the mesonic fields are assumed to be uniform (i.e., space-time 
independent with no quantum fiuctuations) . For the isoscalar vector field, then 



where oOo is space-time independent but depends on pB and 5^ is the Kronecker delta. 
The equation of motion for the vector field specifies cUg '■ 



The equation of motion for a is written for convenience in terms oi y = x/xg, and is of 
the form 



(3.3) 




(3.4) 




(3.5) 
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where M* = yM is the effective mass of the nucleoli and 

Ca = gl/ml; = gl/ml 



(3.6) 



We consider here the normal state of high density matter in which there is no pion 
condensation. 

The diagonal components of the conserved total stress tensor corresponding to the 
Lagrangian (3.1) together with the equation of motion for the fermion field (and a mean 
field approximation for the meson fields) provide the following identification for the total 
energy density (e) and pressure (P) of the many-nucleon system (assumed to be a perfect 
fluid) : 



8c^ 727r4y2 27r2 Jo 



M^{l-y'^Y ^'fcj^^ 7 fkF dkk 



j'^ p,) 

Jo . ,,^on1/2 



8c^ 727rV STT^io ^-^ +M*2)^/2 

The energy per nucleon is 



E = ^" ^" ^^3^^ + + 7^ / dk^ik + M^'^f. (3.8) 



^^c^kl ^ 127r2y2 + k% Jo 

For pure neutron matter 7 = 2 and for nuclear matter 7 = 4. A specification of the 
coupling constants Co-, c^^ now specifies the equation of state. 



3.3 Nuclear matter equation of state 

For nuclear matter, we fix and by fits to two nuclear matter properties : the 
saturation density (n^) and the binding energy per particle at pb = rig- For these we 
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choose the values 0.153 fm ^ and —16.3 MeV respectively, as suggested from analysis of 
experimental data (Moller, Myers, Swiatecki and Treiner This gives 

c„ = 6.2033 fm^ 
= 2.9378 fm^. (3.9) 

This leads to a value of 0.78M for the effective mass of the nucleon in saturating nuclear 
matter. The value of K at saturation density that we get is ~ 700 MeV. 
Values of the various thermodynamical quantities for nuclear matter and neutron matter 
for various densities are presented in Table 3.1a and Table 3.1b. The quantity fi appearing 
in this table is the chemical potential, given by the relationship : = {P + e) / ps and p 
is the total mass-energy density, e/c^. 

In heavy-ion collision experiments, hot hadronic matter is produced at temperatures 
upto 100 MeV, which contain upto 25% of their energy in nuclear resonances and mesonic 
degrees of freedom. After making allowance for (model-dependent) corrections for the 
thermal part of the energy, first estimates of the energy per nucleon of nuclear matter 
(for ksT = 0) have been made (Stock [^] for a discussion). These estimates by Stock 



| 57| are unlikely to be firm, however, for our purpose we have considered here as guide 
for comparison with our results. In Fig. 3.1, we present a comparison of available such 
estimates with the prediction of the equation of state considered by us here. Upto a 
density of iug, there is satisfactory agreement between the two. 
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TABLE 3.1a 

EQUATION OF STATE OF DEGENERATE NUCLEAR MATTER 
AS GIVEN BY THE PRESENT MODEL 



kp Pb y P P E 11 



{f'>T^ ^) ^ 9 cm ^ dyn cm ^ MeV MeV 
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0. 
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.2 
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,92 
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.79 


2. 


,44 


E 


14 
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972.47 


1125.00 


1 


.8 


0. 


,394 


0, 


.74 


7. 


,04 


E 


14 


1.23 E 35 


1002.30 


1197.90 


1 


.9 


0. 


,463 


0, 


.76 


8. 


,57 


E 


14 


1.75 E 35 


1037.30 


1272.40 


2, 


.0 


0. 


,540 


0, 


.78 


1. 


,04 


E 


15 


2.35 E 35 


1076.30 


1347.60 


2, 


.1 


0. 


,625 


0, 


.80 


1, 


,25 


E 


15 


3.05 E 35 


1118.40 


1423.00 


2, 


.2 


0. 


,719 


0, 


.83 


1, 


,49 


E 


15 


3.86 E 35 


1163.10 


1498.40 


2, 


.3 


0. 


,822 


0, 


.86 


1. 


,77 


E 


15 


4.79 E 35 


1209.70 


1573.80 


2, 


.4 


0. 


,934 


0, 


.88 


2. 


,09 


E 


15 


5.85 E 35 


1257.90 


1649.00 


2, 


.5 


1. 


.055 


0, 


.91 


2, 


,16 


E 


15 


7.05 E 35 


1307.10 


1721.00 



37 



TABLE 3.1b 



EQUATION OF STATE OF DEGENERATE NEUTRON MATTER 
AS GIVEN BY THE PRESENT MODEL 



k 


F 


PB 




y 











P 


E 


/i 










g cm 


-3 


dyn cm'"^ 


MeV 


MeV 


1 


.0 


0.034 


0, 


.95 


0, 


.57 


E 


14 


-9.19 E 31 


941.03 


939.33 


1 


.1 


0.045 


0, 


.94 


0, 


.75 


E 


14 


-2.05 E 32 


940.39 


937.53 


1 


.2 


0.058 


0, 


.92 


0, 


.98 


E 


14 


-3.76 E 32 


939.49 


935.47 


1 


.3 


0.074 


0, 


.90 


1, 


.24 


E 


14 


-5.78 E 32 


938.41 


933.54 


1 


.4 


0.093 


0, 


.87 


1, 


.55 


E 


14 


-6.91 E 32 


937.32 


932.67 


1 


.5 


0.110 


0, 


.84 


1, 


.90 


E 


14 


-3.56 E 32 


936.57 


934.62 


1 


.6 


0.140 


0, 


.80 


2, 


.31 


E 


14 


1.26 E 33 


936.82 


942.49 


1 


.7 


0.170 


0, 


.77 


2, 


.78 


E 


14 


5.66 E 33 


939.12 


960.40 


1 


.8 


0.200 


0, 


.75 


3, 


.32 


E 


14 


1.46 E 34 


944.77 


991.16 


1 


.9 


0.230 


0, 


.73 


3 


.94 


E 


14 


2.93 E 34 


954.84 


1033.80 


2, 


.0 


0.270 


0, 


.73 


4, 


.67 


E 


14 


4.99 E 34 


969.73 


1084.90 


2, 


.1 


0.310 


0, 


.73 


5, 


.52 


E 


14 


7.60 E 34 


989.25 


1141.00 


2, 


.2 


0.360 


0, 


.74 


6, 


.49 


E 


14 


1.08 E 35 


1012.90 


1199.90 


2, 


.3 


0.410 


0, 


.75 


7, 


.62 


E 


14 


1.45 E 35 


1040.10 


1260.30 


2, 


.4 


0.470 


0, 


.77 


8, 


.91 


E 


14 


1.88 E 35 


1070.20 


1321.50 


2, 


.5 


0.530 


0, 


.79 


1, 


.04 


E 


15 


2.37 E 35 


1102.70 


1383.20 



NOTE.- The respective columns stand for Fermi momentum, nucleon number density, the 
nucleon effective mass factor, the total mass-energy, the pressure, the energy per nucleon 
and the nucleon chemical potential. The numbers following the letter E represent powers 
of ten in all the tables. 
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Fig. 3.1: First generation of estimates, from heavy-ion collision data, of energy per 
nucleon, E, of nuclear matter plotted against ub = Pb (in units of Ug). The crosses corre- 
spond to pion data and the circles to radial energies (Stock for a detailed discussion). 
The dashed curve corresponds to the present model (nuclear matter). 



39 



3.4 Neutron star matter equation of state 

At high densities typical of interiors of neutron stars, the neutron chemical potentials 
will exceed the combined masses of proton and electron. Asymmetric nuclear matter 
with an admixture of electrons (rather than pure neutron matter) is, therefore, a more 
likely composition of matter in neutron star interiors. The concentrations of protons and 
electrons (denoted by Up and rif. respectively) can be determined using conditions of beta 
equilibrium {n ^ p + e + i?) and electrical charge neutrality : 

Up = He (3.10) 

{Hi = chemical potential of particle species i). 

Since nuclear force is known to favour isospin symmetry, and since the symmetry 
energy arising solely from the Fermi energy is known to be inadequate to account for the 
empirical value of the symmetry energy (~ 32 MeV), we include the interaction due to 
isospin triplet p- meson in Eq. (3.1) for purpose of describing neutron-rich matter. That 
is, we add the following terms : 

to the right hand side of Eq. (3.1) in order to describe the asymmetric matter. Here, ~p ^ 
stands for the p- meson field with mass mp, Qp is the couphng strength and 

G^u = d^~p^ - du~p^. (3.12) 

Strictly speaking, the p- meson should couple to the total conserved current (Glendenning, 
Banerjee and Gyulassy l6l|]). In the above, we have coupled the p- meson to the baryons, 
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which are not the only possible source of isospin current. However, for the ground state 
equation of state, in the mean-field approximation, only the baryon part of the isospin 
current will survive (Glendenning |}41|). 

The equation of motion for ~p in the mean field approximation where ~p ^ is replaced by 
its uniform value (here superscript 3 stands for the third component in isospin space) , 
gives 

T^yo = \9p E <^7o^^^>B, (3.13) 

B=n,p 

where the sum is over neutrons and protons. This gives the following density dependence 
for the field variable : 



Po=|^K-^n)- (3.14) 



The symmetric energy coefficient that follows from the semi-empirical nuclear mass for- 
mula (that is, the coefficient of the term {np — nn)'^/{np + n„)^ in the mass formula). 



is: 



where Cp = gp/rUp and kp = {Gtt'^Pb/iY^^ {pb = np + Un) ■ We fix the coupling constant 
Cp by requiring that agym correspond to the empirical value 32 MeV. This gives 

Cp = 4.6617 /ml (3.16) 

It is noted that the p- meson will contribute a term = m'^{pl)^/2 to the energy density 
and pressure. Table 3.2 lists the pressure versus the total mass-energy density for the 
neutron-rich matter in beta equilibrium. 
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TABLE 3.2 



PRESSURE VS DENSITY FOR NEUTRON-RICH MATTER 



P P 



{g cm ^) {dyn cm ^) 



4. 


,531 


E 


15 


1, 


.619 


E 


36 


3, 


,686 


E 


15 


1, 


.285 


E 


36 


2. 


,695 


E 


15 


8, 


.978 


E 


35 


2. 


,109 


E 


15 


6, 


.719 


E 


35 


1. 


,784 


E 


15 


5, 


.472 


E 


35 


1. 


,504 


E 


15 


4, 


.406 


E 


35 


1. 


,211 


E 


15 


3 


.296 


E 


35 


1. 


,017 


E 


15 


2, 


.565 


E 


35 


9. 


,116 


E 


14 


2, 


.171 


E 


35 


8. 


,175 


E 


14 


1 


.822 


E 


35 


7. 


,028 


E 


14 


1, 


.400 


E 


35 


6. 


,059 


E 


14 


1, 


.051 


E 


35 


5. 


,576 


E 


14 


8, 


.813 


E 


34 


5. 


,037 


E 


14 


6, 


.967 


E 


34 


4. 


,561 


E 


14 


5, 


.410 


E 


34 


4, 


,063 


E 


14 


3 


.887 


E 


34 


3. 


,564 


E 


14 


2, 


.534 


E 


34 


3. 


,018 


E 


14 


1, 


.352 


E 


34 
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3.5 Results 



Using the chiral sigma model and adopting the approach that the isoscalar vec- 
tor field, needed to provide saturating binding energy of degenerate nuclear matter be 
generated dynamically, we have obtained an equation of state of degenerate nuclear and 
neutron-rich matter at high densities. 

The maximum mass of the neutron stars is calculated by integrating the structure 
equations, which will be discussed in the chapter 7 in detail. A comparison with the 
maximum mass of (nonrotating) neutron stars predicted by others, using models based on 
recent field theoretical equations of state is given in Table 3.3. The choice of the equation- 
of-state models in Table 3.3 is representative, but by no means exhaustive. Included 
in this comparison are equation of state models due to Alonso and Cabanell ||3^ and 



Prakash and Ainsworth E3] which are also based on the sigma model, but differ from 



comes 



our model in the details. The equation of state (II) of Alonso and Cabanell pTj . 
from a determination of the free parameters of the linear sigma model with an explicit 
symmetry-breaking term and is coupled to uo- and p- mesons, in a renormalizable way. 



Prakash and Ainsworth included the sigma meson one-loop contributions, but the 
isoscalar vector field was not generated dynamically, so that its role is reduced to an 
empirical one, allowing for arbitrary variations in its coupling constant (thereby making 
it possible to obtain any desired value of the nuclear matter compression modulus). The 
vector field plays no role in determining the value of the effective mass of the nucleon in 
such an approach. A comparison of this equation of state with our model for neutron-rich 
matter and pure neutron matter is shown in Fig. 3.2. The equation of state of Serot 

43 



TABLE 3.3 



MAXIMUM MASS NEUTRON STARS FOR FIELD THEORETICAL 
EQUATION OF STATE MODELS 



Equation of state Ref 


erence 


Us 


K 




R 












MeV 




km 


Serot 


'62 






0.193 


540 


2.54 


12.28 


Chin BSI 






0.193 


471 


2.10 


10.57 


Alonso-Cabanell [B^ equation of state II 


0.172 


225 


1.94 


10.90 


Prakash-Ainsworth 
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isl = 16.27) 


0.160 


225 


1.83 


9.89 


Glendenning 






0.153 


300 


1.79 


11.18 


Present model fc5 


0.153 


700 


2.59 


14.03 



NOTE.- n^ and K are nuclear matter saturation density and compression modulus. 



62| is calculated in the mean- field approximation in the Walecka model including cr— , uj- 



and p- mesons and that of Chin p3[ with one-loop corrections in the a — uj model. The 



equation of state due to Glendenning |^1[ is along similar lines, and includes one-loop 
corrections and also scalar self-interactions (upto quartic order), whose magnitudes are 
adjusted to reproduce empirical saturation properties. This work takes into account the 
effect of hyperons in beta equilibrium in addition to electrons and muons. 

Table 3.3 implies that the present neutron-rich matter equation of state is compara- 
tively 'stiff' as far as neutron stars are concerned. This is reflected in the value of the 
maximum mass of neutron stars, which is the largest for the present model. It may be 
mentioned here that observational evidence in favour of a stiff equation of state comes 



from the identiflcation by Triimper et al. of the 35 day cycle of the pulsating X-ray 



source Her X-1 as originating in free precession of the rotating neutron star (Pines 

In constructing the neutron star matter equation of state, we have restricted ourselves 
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Fig. 3.2: Pressure versus total mass-energy density. Curve 1 corresponds to neutron- 
ricli matter in beta equilibrium (including the contribution from p— meson exchange) : 
present model. Curve 2 corresponds to pure neutron matter : present model. Curve 3 
(dashed curve) is due to Prakash and Ainsworth P2 . 
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to (n, p, e~) matter in beta equilibrium. At high densities, hyperons can also appear 
(Glendenning Kapusta and Olive |3^); this is expected to reduce the 'stiffness' in 
the equation of state, and a consequent reduction in the maximum mass of neutron stars. 
Another point that we have not investigated here is the possible role of the p— meson 
tensor interaction as far as the symmetry energy is concerned. 
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Chapter 4 

Astrophysical application I : 
Structure and radial oscillation of 
nonrotating neutron stars 

4.1 The Neutron Star Structure 

The structure of a neutron star is characterized by its mass and radius. Additional 
parameters of interest are the moment of inertia and the crust thickness. These are 
important for the dynamics and transport properties of pulsars. 

The space-time for a spherically symmetric gravitating system is described by the 
Schwarzschild metric 

ds^ = e^c^de - r\d9^ + sm^ddcj?) - e^dr^, (4.1) 



where z/, A are functions of r only [117]. Corresponding to this space-time metric, the 
equations that describe the hydrostatic equilibrium of degenerate stars without rotation 
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in general relativity {i.e., where temperature and convection are not important) can be 
written as 11 1711 : 



dp _ G{p + p/c^){m + Attt^p/c^) ^^^^ 
dr r2(l — 2Gm/rc^) 



dm 9 

— = 47rrV 4.3 
dr 

where p and p are the pressure and total mass energy density and m(r) is the mass 
contained in a volume of radius r. Given an equation of state p(p), Eqs. (7.2) and (7.3), 
called Tolman-Oppenheimer-Volkoff (TOV) equation, can be numerically integrated, for 
a given central density, to obtain the radius R and gravitational mass Mq = m{R) of the 
star. 

The moment of inertia I of the rotating neutron star, rotating with angular velocity 
Q as seen by a distant observer is given by [p.l8| 



/ = l^f^) (44) 

where a;(r), the angular velocity of the star fluid relative to the local inertial frame, is 
given by 

d / . duj\ ■} dj , ^ 

-^ir'3-r)+^r^^f = Q, 4.5 
dr^ dr ^ dr 

where 

j(r) = 6-^^/2(1 - 2GMG/r^f'^ (4.6) 
and satisfies the boundary conditions 
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(^)„« = 0;-(oo) = a (4.7) 

Here the potential function i^(r), relating the element of proper time to the element of 
time at r = oo is given by 

du _ 2G (m + Anr^P/c^) 

d^ ~ r^d^ (1 - 2Gm/rc^) ' ^'^ 

The definition (7.4) for I includes the relativistic Lens-Thirring effect to order Q^. For 
neutron star, the relativistic effect are important, and so one gets / as defined by Eq. 
(7.4) to be in excess of classical definition (2/5) Mi?^. 

For the numerical integration to obtain the structure parameters, it is sufficient to start 
with an arbitrary value of iy{0), which is then rescaled to satisfy the surface condition 

.(fl)=Kl-^). (4.9) 

SO that i/(oo) = 0. Likewise, a;(0) is initially chosen to be an arbitrary constant, and a 
value of Q given by 

"=-(fi)4«(f)~« ("0) 

obtained. A new starting value LJnewi^) corresponding to any specified flnew is given by 

(0) — nnew<^(0)/n. (4-11) 

To integrate the TOV equations, one needs to know the equation of state p(p), for 
the entire expected density range of neutron star, starting from the high density at the 
center to the surface densities. The composite equation of state for the entire neutron star 
density span, was constructed by joining the equation of state of high density neutron-rich 
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TABLE 7.1 



PRESSURE VS DENSITY FOR NEGELE AND VAUTHERIN 
EQUATION OF STATE 



P P 



{g cm ^) {dyn cm ^) 



1 


586 


E 


14 


8 


617 


E 


32 


9 


826 


E 


13 


3 


807 


E 


32 


6 


193 


E 


13 


1 


835 


E 


32 


3 


767 


E 


13 


8 


564 


E 


31 


2 


210 


E 


13 


3 


789 


E 


31 


1 


496 


E 


13 


2 


095 


E 


31 


9 


611 


E 


12 


1 


095 


E 


31 


6 


248 


E 


12 


6 


184 


E 


30 


3 


833 


E 


12 


3 


621 


E 


30 


2 


202 


E 


12 


2 


276 


E 


30 


1 


471 


E 


12 


1 


694 


E 


30 


9 


728 


E 


11 


1 


228 


E 


30 


6 


610 


E 


11 


8 


633 


E 


29 


5 


228 


E 


11 


6 


741 


E 


29 



matter (that we have discussed in chapter 3) to that given by (a) Negele and Vautherin 
| 119|| for the density region (10^^ — 5 x 10^°) g cm~^, (b) Baym, Pethick and Sutherland 



| 120|| for the region (5 x 10^° — 10^) g cm^"^ and (c) Feynman, Metropolis and Teller [ |121 
for p < 10^ g cm~^. These densities regions are tabulated in Table 7.1-7.3. 

For a given equation of state -P(p), and a given central density p(r = 0) = pc, the Eqs. 
(7.2-7.3) are integrated numerically with the boundary condition : 

m{r = 0) = (4.12) 

to give R and Mq. The radius R is defined by the point where P ~ 0, or, equivalently, 
p = ps, where ps is the density expected at the neutron star surface (about 7.8 g cm~^). 
The total gravitational mass is then given by : Mq = m{R). 
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TABLE 7.2 



PRESSURE VS DENSITY FOR BAYM, PETHICK AND SUTHERLAND 

EQUATION OF STATE 



P P 



{g cm ^) {dyn cm ^) 



5, 


.254 


E 


10 


5. 


.949 


E 


28 


3, 


.313 


E 


10 


3. 


,404 


E 


28 


2, 


.090 


E 


10 


1. 


,938 


E 


28 


1, 


.318 


E 


10 


1, 


,048 


E 


28 


8, 


.312 


E 


09 


5. 


,662 


E 


27 


4, 


.164 


E 


09 


2. 


,356 


E 


27 


1, 


.657 


E 


09 


6. 


,860 


E 


26 


1, 


.045 


E 


09 


4. 


,129 


E 


26 


5, 


.237 


E 


08 


1. 


,629 


E 


26 


2, 


.624 


E 


08 


6. 


,676 


E 


25 


6, 


.589 


E 


07 


1. 


,006 


E 


25 


1. 


.655 


E 


07 


1. 


,435 


E 


24 


6, 


.588 


E 


06 


3. 


,911 


E 


23 


1, 


.044 


E 


06 


2. 


,318 


E 


22 


1, 


.654 


E 


05 


1. 


,151 


E 


21 


2, 


.622 


E 


04 


4. 


,968 


E 


19 


1 


.044 


E 


04 


9. 


,744 


E 


18 



TABLE 7.3 

PRESSURE VS DENSITY FOR FEYNMAN, METROPOLIS AND TELLER 

EQUATION OF STATE 



P 

{g cm~^) 
2.120 E 03 
4.500 E 01 
1.640 E 01 
1.160 E 01 
8.150 E 00 
7.900 E 00 



P 

{dyn cm^^) 
5.820 E 15 
1.700 E 14 
1.400 E 13 
1.210 E 12 
1.010 E 11 
1.010 E 10 
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We solved the TOV equation numerically using predictor-corrector method, which 
gives better accuracy compared to existing results of the gravitational mass and the radius. 
In this method, we use logarithm of pressure and energy density e.g., equation of state to 
take care of lower energy density and pressure, because at low densities the energy density 
falls very rapidly at the surface of the star. 

Table 7.4 lists neutron star structure parameters as predicted by our equation of state 
(chiral sigma model) for neutron-rich matter, which has already been discussed in chapter 
3. The maximum gravitational mass for stable non- rotating neutron star predicted by 
our model is 2.59 Mq (Mq = solar mass). This occurs for a central density of 1.4 x 10^^ 
g cm~^. The corresponding radius and entire crustal length A of the star are 14.03 km 
and 1 km respectively. The crust length A is defined as the distance overwhich the density 
falls from p = 2.4 x 10^^ to 7.8 g cm~^ i.e., the surface. Fig. 7.1 and Fig. 7.2 show the 
plots of mass vs central density and moment of inertia vs mass. The maximum moment 
of inertia is 4.79 x 10^^ g cm^. 

Observationally, masses of neutron stars are estimated from compact binary systems, 
one member of which is a pulsar. The most precise estimate comes from the pulsar PSR 
1913+16, which gives (1.442 ±0.003) Mq. A recent compilation of the estimated masses by 
X-ray pulsars (Nagase [p.22|| ) gives the maximum mass (corresponding to Vela X-1 pulsar) 
to be (1.77 ±0.21) Mq. Stable neutron star masses predicted by the present equation of 
state discussed in chapter 3 are thus compatible with the observational estimates. The 
surface red shift factor provides a probe for neutron star structure, if one presumes that 
observed 7— ray bursts are gravitationally red shifted e^e~ annihilation lines, produced 
near their surface. The surface red shift ratio (a) defined as 
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Fig. 7.1: Gravitational mass (Mq) of non- rotating neutron stars versus central density 
(pc) as predicted by the present model equation of state (neutron-rich matter). The 
maximum stable mass is 2.59 Mq. The corresponding central density, Pc — 1.5 x 10^^ 
g cm~^ and the radius is 14.0 km. 
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Fig. 7.2: Neutron star moment of inertia (/) versus gravitational mass (Ma), as predicted 
by the present model. 
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TABLE 7.4 



NEUTRON STAR STRUCTURE PARAMETERS FOR THE PRESENT MODEL 
NEUTRON STAR EQUATION OF STATE 



Pc 


R 


M/Mq 


/ 


a 


A 




{km) 




{g cm^) 




{km) 


3.5 E 14 


14.30 


1.02 


1.31 E 45 


0.889 


4.06 


4.0 E 14 


14.73 


1.36 


2.05 E 45 


0.853 


3.11 


4.5 E 14 


14.99 


1.64 


2.73 E 45 


0.822 


2.55 


5.0 E 14 


15.13 


1.86 


3.30 E 45 


0.798 


2.18 


6.0 E 14 


15.18 


2.16 


4.07 E 45 


0.761 


1.76 


8.0 E 14 


14.94 


2.45 


4.70 E 45 


0.719 


1.36 


1.0 E 15 


14.62 


2.55 


4.79 E 45 


0.696 


1.18 


1.2 E 15 


14.30 


2.59 


4.70 E 45 


0.682 


1.07 


1.4 E 15 


14.03 


2.59 


4.54 E 45 


0.674 


1.00 



NOTE.- The respective columns stand for central density, radius, gravitational mass, 
moment of inertia calculated for angular velocity = {GMg/ R^)^^'^, and the surface red 
shift ratio a. 

a = {l-2GMG/RcY^^ (4.13) 
is expected to be (0.78 ±0.02) on the basis of observed data (Friedman and Pandharipande 



26| for a discussion). The present neutron-rich matter equation of state, which has already 
been discussed, gives for a 1.4 Mq neutron star : R = 14.77 km, I = 2.15 x lO'*^ g cm^ 
A=3.0 km and the red shift ratio (at the surface) = 0.85. The corresponding central 
density is 4.06 x 10^^ g cm~^. 

4.2 Neutron stars radial oscillations 



Since the original suggestion by Cameron [ |123[ that vibration of neutron star could 



excite motions, which might have interesting astrophysical applications, there has been 
a series of investigations into the vibrational properties of neutron stars. The earliest 
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detailed calculations, carried out by Meltzer and Thorne ||124|| and Thorne ||125|| examined 
the radial as well as nonradial oscillations of neutron stars, using then available equation 
of state, such as the Harrison- Wakono- Wheeler |126| equation of state. These and other 



early studies, e.g., Wheeler [ 127 |; Chau | 128 | and Occhionero [129], indicated that the 
fundamental mode radial oscillation periods for the neutron stars would typically lie in 
the vicinity of about 0.4 ms, and that the first few quadrupole oscillations would have 
periods that are also fractions of a millisecond. Furthermore, these oscillation periods 
were estimated to be damped by gravitational radiation with damping time scales of the 
order of one second. Also, there are recent papers on nonrotating, rapidly rotating neutron 
stars and slowly rotating stars. For example, in 1990, Cutler, Lindblom and Splinter [|130|| 
computed the frequencies and damping times due to viscosity and gravitational radiation 
for the lowest frequency modes of a wide range of nonrotating fully relativistic neutron 



star models. In a subsequent paper. Cutler and Lindblom ||131|| reviewed and extended 
the formalism for computing the oscillation frequency of rapidly rotating neutron stars in 



the post-Newtonian approximation. Recently, Kojima |p.32|] calculated the rotational shift 
of normal frequencies in polytropic stellar models in the framework of general relativity. 
The stellar rotation is assumed to be slow and first-order rotational effects are included to 
the eigenfrequencies of the nonrotating stars. In this section, we present the calculation of 
radial oscillation periods of neutron stars, using our equation of state based on the chiral 
sigma model. 

The equation governing infinitesimal radial pulsations of a nonrotating star in general 



relativity was given by Chandrasekhar ||133|| , and it has the following form : 

ar'^ ar 



(4.14) 
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where ^(r) is the Lagrangian fluid displacement and ccr is the characteristic eigenfrequency 
(c is the velocity of light). The quantities F, G, H depend on the equilibrium profiles of 
the pressure and density of the star, and are given by 



F = -e"^eTp/(p + pc') (4.15) 



Pf + r|}/(p + pc^) (4.16) 



p + pc^ ar p + pc^ ^ 
where F is the adiabatic index, defined in the general relativistic case as 



and 



^_d\Tp ^2pdT ^ 2T dp _ 2Tp 
dr r r dr r dr 

l^(^rp + 2p— + 2r— - —) 

A dr^dr dr dr r ' 

1^ fdi'\2 1 d^v 



-dr^ 2 dr'^ 

The boundary conditions to solve the pulsation equation (7.14) are 



e(r = 0) = (4.20) 



dn p^/"^ r) 

Sp{r ^R) = -i£ - Tp-^Q-ye-^''i)\r^n - (4.21) 
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Here T is the adiabatic index. Since p vanishes aX, r = R, it is generally sufficient to 
demand 



Equation (7.14), subject to the boundary conditions Eq. (7.20) and (7.21), is a Sturm - 
Liouville eigen value equation for cr^. 

The following results follow from the theory of such equations : 

1. The eigen values are all real. 

2. The eigen values from an infinite discrete sequence 



3. The eigen function corresponding to has no nodes (fundamental mode) in the 
interval < r < R; more generally, has nodes in this interval. 

4. The are orthogonal with weight function cor'^ : 



5. The $,n form a complete set for the expansion of any function satisfying the boundary 
condition Eqs. (7.20) and (7.22). 

An important consequence of item (2) is the following : 

If the fundamental radial mode of a star is stable (cr^ > 0), then all the radial modes are 
stable. Conversely, if the star is radially unstable, the fastest growing instability will be 
via the fundamental mode (cr^ more negative than all other cr^). 



^ finite at r — R. 



(4.22) 



a„ < (Ti < . . . < (j„ < 



(4.23) 
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We solved [ p,34|| Eq. (7.14) for the eigenvalue a by writing the differential equation 
as a set of difference equations. The equations were cast in tridiagonal form and the 
eigenvalues were found by using the EISPACK routine. This routine finds the eigenvalues 
of a symmetric tridiagonal matrix by the implicit QL method. 

Results for the oscillations of neutron star corresponding to chiral sigma equation of 
state are illustrated in Fig. 7.3. The plot in Fig. 7.3 are for the oscillation time period 
(= 27r/ccr) versus the gravitational mass Mq- The fundamental mode and the first four 
harmonics are considered. The period is an increasing function of Mg, the rate of increase 
being progressively less for higher oscillation modes. The fundamental mode oscillation 
periods for neutron star are found to have the following range of values : (0.35 - 0.50) 
milliseconds. For higher modes, the periods are < 0.2 milliseconds. We shall compare 
our results with those for a different choice of the equation of state, namely, that given by 



Wiringa et al. [|lT|. We also perform a similar calculation for strange quark stars, using a 



realistic equations of state. 



4.3 Radial oscillations of quark stars 



There are strong reason for believing that the hadrons are composed of quarks, and 
the idea of quark stars has already existed for about twenty years | ]135| , |136| |. Calculations 
of the possible phase transition from baryon matter to quark matter in models of cold, 
compact stars have been performed by several groups (|9^, |137| , P^ , |138| |), but the results 
are not conclusive concerning the existence of quark matter inside neutron stars. In 1984, 
it was suggested that strange matter i.e., quark matter with strangeness per baryon of 
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Fig. 7.3: Periods of radial pulsations as functions of the gravitational mass for our 
equation of state (chiral sigma model). The labels 1, 2, 3, 4, 5 correspond respectively to 
the fundamental and the first four harmonics. 
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order unity, may be the true ground state [^. The properties of strange matter at zero 
pressure were subsequently examined, and it was found that strange matter can indeed be 
stable for a wide range of parameters in the strong interaction calculations P^. Details 



of the extension to finite pressure and the so-called strange stars are given in |^8|, |139 
The problem of the existence of strange stars is, however, still unresolved ( [|100|| ). 

An important question is how can one possibly distinguish between quark stars and 
neutron stars. It has been suggested to use measurements of the surface gravitational red 
shift z (Schwarzschild), since different equation of state gives different results for z{Mg\ 



M being the star mass ||140| , |141|| . The region of allowed high density equation of state 
may be narrowed further by the observations of pulsar periods [ p.42| , |143|| . Given a sub- 
millisecond pulsar, we may argue that the ability of such a fast rotating star to avoid 
rotational break-up induces severe restrictions and a conventional neutron star will not 
be able to resist the large centrifugal forces. The problem of rapid rotation of compact 
stars receives much attention [|144|, |145| , |146| , |14711 , and although no sub-millisecond pulsar 
is seen among the about 500 pulsars observed so far, further observations may well reveal 
such an object. Also, among the criteria suggested for distinguishing quark star from 
neutron star are the neutrino cooling rate [|14^ , p.49| , |150|| , transport properties such as 
bulk viscosity [|151| , |152|| , and sub millisecond period rotation rates ||145|] . 



These so-called strange stars have rather different mass-radius relationship [ [L39| | than 
neutron stars, but for stars of mass = 1.4 M©, the structure parameters of quark stars are 
very similar to those of neutron stars. Since pulsars are believed to be (rotating) neutron 
stars, and since available binary pulsar data suggest their masses to be close to 1.4 M©, 
it has been conjectured that at least some pulsars could be quark stars. 
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Recently, Haensel et al. ||153|| have emphasized that pulsation properties of a neutron 
star can yield information about the interior composition, namely, whether the interior 
has undergone a phase transition to quarks. The main idea is to know the damping times. 



which will be modified if there is a quark matter core. In their study, Haensel et al. [153] 
used polytropic model for the equation of state for nuclear matter as well as quark matter, 
and the Newtonian pulsation equation to calculate the eigenfrequencies. The strange 



quark mass and the quark interactions are important for the structure of quark stars ||99 
This suggests that the equation of state of strange quark matter will have a role to play in 
determining the pulsation features of quark stars. Clearly, for a more exact understanding 
of the vibrational properties of quark stars, use of realistic equation of state for quark 



matter, and the general relativistic pulsation equation, are desirable. Cutler et al. [130| 
have calculated the frequencies and damping times of radial pulsations of some quark star 
configurations, using the general relativistic pulsation equation, but for quark matter, they 
adopted the MIT bag model in its simplest form, namely, non-interacting and massless 



quarks. The purpose of this work [|134|| is to calculate the range of eigenfrequencies of 
radial pulsations of stable quark stars (using the general relativistic pulsation equation) 
and to investigate the sensitivity of the eigenfrequencies on the equation of state. 

The equation of state used by us incorporates short-range quark-gluon interactions 
perturbatively to second order in the coupling constant ac- The long-range interactions are 
taken into account phenomenologically by the bag pressure term B. We incorporate the 
density dependence of by solving the Gell-Mann-Low equation for the screened charge. 
The parameters involved are the strange quark mass m^, B and, the renormalization point 
yUo, obtained by demanding that the bulk strange matter be stable at zero temperature 
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and pressure, with energy per baryon less than the lowest energy per baryon found in 
nuclear matter. For completeness, we also do the calculations for the MIT bag model. 

At high baryonic densities, bulk strange matter is in an overall colour singlet state, 
and can be treated as a relativistic Fermi gas interacting perturbatively. The quark 
confinement property is being simulated by the phenomenological bag model constant B. 
Chemical equilibrium between the three quark fiavours and electrical charge neutrality 
allow us to calculate the equation of state from the thermodynamic potential of the system 
as a function of the quark masses, the bag pressure term B and the renormalization point 
/ic,. To second order in Oc, and assuming u and d quarks to be massless, the thermodynamic 
potential is given by ||154|| : 



n = nu + nd + ns + a„t. + n^, (4.24) 

where fli {i = u, d, s, e) represents the contributions of u, d, s quarks and electrons 
and Qint is the contribution due to interference between u and d quarks and is of order 

: Expressions for Qi and Qint are already given in chapter 5. 
The total energy density and the external pressure of the system are given by 

e = n + B + Y.l^i^i (4-25) 

i 

and 

p=-n-B, (4.26) 

where rii is the number density of the i-th particle species. For specific choices of the 
parameters of the theory (namely, m^, B and /Iq), the equation of state is now obtained 
by calculating e and p for a given value of /i : 
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fl = jj^d = jJ's = jJ'U + /ie, (4.27) 

by solving for /Xe from the condition that the total electric charge of the system is zero. 

There is an unphysical dependence of the equation of state on the renormalization 
point fio, which, in principle, should not affect the calculations of physical observables if 



the calculations are performed to all orders in Uc ||139| , |103|| . In practice, the calculations 
are done perturbatively and, therefore, in order to minimize the dependence on /io the 
renormalization point should be chosen to be close to the natural energy scale, which 
could be either /i^ ~ B^/^ or the average kinetic energy of quarks in the bag, in which 
case, fjjQ ~ 313 MeV. In the present study, our choice of /io is dictated by the requirement 
that stable strange matter occurs at zero temperature and pressure with a positive baryon 
electric charge |154||. This leads to the following representative choice of the parameter 



values: 

Equation of state model 1 : B = 56 MeV fm-^; m, = 150 MeV; /i^ = 150 MeV. 
Equation of state model 2 : B = 67 MeV fm'^' m^ = 150 MeV; = 0. 
Model 2 corresponds to no quark interactions, but a non-zero mass for the strange quark. 
In the limit, m^ — > and etc — > 0, the equation of state has the analytical form 

p=l(e-4i?), (4.28) 

where e is the total energy density. Eq. (7.28) is the MIT bag model. It is independent 
of the number of quark flavours. 

Numerical values of pressure p and total mass-energy density p = e/c^ for the quark 
matter equation of state models used here are listed in Table 7.5. For the sake of compari- 
son, we have included in this table, the equation of state corresponding to non-interacting, 
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TABLE 7.5 

EQUATIONS OF STATE FOR DEGENERATE STRANGE QUARK MATTER 



p(10^^ g cm ^) 


P{10^^ dyn cm~^) 
model 1 


model 2 


MIT bag {B = 56 MeV fm-^) 


6.0 


4.44 


2.23 


6.01 


8.0 


10.13 


7.83 


12.00 


10.0 


15.88 


13.49 


17.99 


12.0 


21.63 


19.17 


23.99 


14.0 


27.41 


24.88 


29.98 


16.0 


33.20 


30.16 


35.97 


18.0 


39.00 


36.36 


41.96 


20.0 


44.82 


42.12 


47.95 


22.0 


50.64 


47.89 


53.95 


24.0 


56.47 


53.67 


59.94 


26.0 


62.30 


59.46 


65.93 


28.0 


68.14 


65.26 


71.92 


30.0 


73.98 


71.06 


77.91 


32.0 


79.83 


76.87 


83.90 


36.0 


91.53 


88.51 


95.89 


40.0 


100.32 


100.16 


107.87 


50.0 


132.58 


129.36 


137.83 



massless quarks as given by the simple MIT bag model with B = 56 MeV fm^^. Among 
these equation of state, the bag model is 'stiffest' followed by models 1 and 2. Equilibrium 
configurations of strange quark stars, corresponding to the above equation of state, are 
presented in Table 7.6, which lists the gravitational mass Mq, radius R, the surface red 
shift z, given by 

z={l- 2GM/c^R)-^/^ - 1 (4.29) 



and the period corresponding to fundamental frequency VLo defined as 130 1 : 



fio = (3G'M/4i?3)i/2 
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(4.30) 



TABLE 7.6 

EQUILIBRIUM STRANGE QUARK STAR MODELS 



Equation of state 


Pc(10^^ g cm ^) 


M/Mq 


R{km) 


Surface red shift {z) 


Po(ms) 


Model 1 


24.0 


1.958 


10.55 


0.487 


0.488 




20.0 


1.967 


10.78 


0.472 


0.503 




16.0 


1.951 


11.02 


0.448 


0.522 




12.0 


1.864 


11.22 


0.401 


0.548 




8.0 


1.521 


11.02 


0.299 


0.591 




6.0 


0.997 


9.93 


0.192 


0.624 




5.0 


0.485 


7.99 


0.104 


0.646 


Model 2 


24.0 


1.863 


10.09 


0.483 


0.468 




20.0 


1.862 


10.29 


0.465 


0.482 




16.0 


1.829 


10.49 


0.435 


0.500 




12.0 


1.710 


10.62 


0.381 


0.527 




8.0 


1.281 


10.14 


0.263 


0.568 




6.0 


0.645 


8.37 


0.138 


0.600 




5.0 


0.092 


1.18 


0.032 


0.622 


MIT Bag 


24.0 


2.021 


10.81 


0.493 


0.500 


model 


20.0 


2.033 


11.04 


0.480 


0.514 


{B = 56 


16.0 


2.023 


11.29 


0.450 


0.533 


MeV fm-3) 


12.0 


1.947 


11.52 


0.410 


0.558 




8.0 


1.635 


11.41 


0.310 


0.604 




6.0 


1.150 


10.52 


0.210 


0.636 




5.0 


0.666 


8.98 


0.130 


0.659 
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as functions of the central density pc of the star. 

We calculated the eigenvalue a by solving the Eq. (7.14) as same way as discussed in the 
previous section. 

Results for the oscillations of quark stars corresponding to equation of state models 1 
and 2 are illustrated in Fig. 7.4. To compare, we have included in Fig. 7.4 the results for 
quark stars corresponding to (a) the simple MIT bag equation of state (non-interacting, 
massless quarks and S = 56 MeV fm~^) and (b) neutron stars corresponding to a recently 



given neutron matter equation of state |TT|. The plots in Fig. 7.4 are for the oscillation 
time period (= l-njco) versus the gravitational mass Mq- The fundamental mode and 
the first four harmonics are considered. The period is an increasing function of M, the 
rate of increase being progressively less for higher oscillation modes. The fundamental 
mode oscillation periods for quark stars are found to have the following range of values: 
MIT bag model : (0.14 - 0.32) milliseconds 
Equation of state model 1 : (0.10 - 0.27) milliseconds 
Equation of state model 2 : (0.06 - 0.30) milliseconds 

For neutron stars (model UV14 + UVII, ref.|]lT|), we find that the range of periods for / 
= mode is (0.25 — 0.4) milliseconds, which is slightly less than the neutron star based 
on chiral sigma model equation of state, presented in previous section. For higher modes, 
the periods are < 0.1 milliseconds, similar to the case of quark stars but less than the 
chiral sigma model values. So, the oscillation periods for neutron star based on chiral 
sigma model equation of state is different from quark stars, which is insignificant. 

Inclusion of strange quark mass and the quark interactions make the equation of state 
a little 'softer' as compared to the simple MIT bag equation of state (Table 7.5). This is 
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Fig. 7.4: Periods of radial pulsations as functions of the gravitational mass. The top 
two and bottom left boxes correspond to strange quark stars. The bottom right box is for 
stable neutron stars corresponding to beta-stable neutron matter, model UV14 + UVII, 



ref. [|rT|. The labels 1, 2, 3, 4, 5 correspond respectively to the fundamental and the first 



four harmonics. 
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reflected in the value of the maximum mass of the strange quark star (Table 7.6). For the 
pulsation of quark stars, this gives, for 1=0 mode eigenfrequencies, values as low as 0.06 
milliseconds. The main conclusion that emerges from our study, therefore, is that use of 
realistic equation of state can be important in deciding the range of eigenfrequencies, at 
least for the fundamental mode of radial pulsation. The results presented here thus form 



an improved first step of calculations on the lines presented by Haensel et al. ||153|| , whose 
numerical conclusions are expected to get altered. 

Since we considered the vibrations of neutron stars and quark stars, it is important to 
study the time scale for damping of the vibrations. Regarding the time scale of damping 
of the vibration, here we don't calculate exactly, but the damping time is approximately 
same as Ref. [|155|| , considered by Madsen, because, the oscillation time is taken to be 
10~'^s, which is typical for the fundamental mode in our case. The discussion by Madsen 



155 was based on rather crude estimates. 



69 



Chapter 5 

Astrophysical application II : 
Cooling of neutron stars with quark 
core 

5.1 Introduction 

When a neutron star is formed in the collapse of a stellar core, it rapidly cools down 
by neutrino radiation. The interior temperature drops to about less than 10^° K within 
minutes and to about IQ^ K within one year. Neutrino emission dominates the subsequent 
cooling of the neutron star, until the interior temperature falls to about 10^ K, with a 
corresponding surface temperature of about 10^ K. Thus the photoemission also begin 
to play an important role. 

The cooling curves (observed temperature as a function of time) depend on a number 
of interesting aspects of the physics of neutron star. It turns out that the equation of state 
as well as the mass of the neutron stars do not influence the cooling in a sensitive manner, 
but the possible existence of a superfluid state of the nucleons plays some role, and the 
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existence of a pion condensate or a quark phase in the central region would have dramatic 
effects. For conventional cooling scenarios, neutrino emission dominates the cooling for 
about 10^ years. 

The theoretical cooling models [ 156 , 157 ] have been refined | 158 , 159 1 in recent years by 



a number of authors. The search for thermal radiation of pulsars has so far in most cases 
has led only to upper bounds for the surface temperature, still interesting comparisons 
between theory and observations can be made. 

The so-called standard model of neutron star cooling is based upon neutrino emission 



from the interior that is dominated by the modified URCA process |160 



{■n,p)+p + e ^ (n,p) + n + z/e, (5.1) 

{n, p) + n ^ {n, p) + p + e~ + z/g. (5.2) 

The direct URCA process 

n p + e~ + Ue, (5.3) 

p + e~ n + Ue, (5.4) 

is not usually considered because it is strongly suppressed in degenerate matter because 
of the requirement of energy and momentum conservation. Since this process has recently 
been revived ||161||, we repeat the simple argument. The fermions n, p, e~ participating 



in the process have energies lying within T of the Fermi surface. By energy conservation, 
the neutrino and antineutrino energies are then also ~ T. But the Fermi momenta of the 
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electrons and protons are small compared to the neutron Fermi momentum and thus the 
processes (8.2) are strongly suppressed by momentum conservation. 

This is no longer a case when an additional neutron, which can absorb energy and 
momentum, takes part in the process, as in (8.1), (8.2). A pion condensate would have 
the same effect as that of a spectator neutron. 

It has recently been argued ||161|| that the proton concentration in a neutron star might 



be so high that the momentum conservation, namely 

P/(P) +P/(e) >P/(^), (5.5) 

might be satisfied. For an n, p, e mixture we have Up = rie and thus the condition 
becomes 

Hn < 8np. (5.6) 

The proton fraction x = Up/n, where n = (un + rip) is the total baryon density, is then 
given by 

x>]-- 11.1% (5.7) 
9 

If the electron chemical potential exceeds the muon rest mass = 105.7 MeV, muon 
will also be present in dense matter, and this will increase the threshold proton concentra- 
tion. If fie ^ frifj^, the threshold proton concentration is ~ 0.148; for smaller values of /ig, 
the threshold concentration lies between | and 0.148. At densities typical of the central 
regions of neutron stars, the calculated proton concentration of matter is very sensitive to 
the choice of physical model, and in reality it might exceed the threshold value as Lattimer 
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et al. [|161|] discuss. Estimates of proton fraction as a function of baryon density for a num- 
ber of different equation of states indicate that, (a) the estimated proton concentrations 
depend sensitively on the assumptions made about the microscopic interactions, which 
are poorly known, and (b) it is quite possible that the proton concentrations are large 
enough to allow the direct Urea process to occur. As calculations by Wiringa, Fiks and 



Fabrocini |jTT[ demonstrate, the form of the three-body interaction, especially its isospin 
dependence has a large influence on the proton fraction. 

Let us now estimate the rate at which antineutrino energy emitted per unit volume in 
the reaction given by (8.3). This may be done by using Fermi's "golden rule". Neglecting 
for the moment the effects of possible superfiuid of neutrons and superconductivity of 
protons, one finds 



^/3 = y2^GW^c(l + 3<7i) 
xr2i(l - n2)(l - ^3) X ei5^{pi - P2 - Ps - Pi), (5. 



where rii is the Fermi function and the subscript i=l to 4 refer to the neutron, proton, 
electron, and antineutrino respectively. The pt are four-momenta, and £4 is the antineu- 
trino energy. The sum over states is to be performed only over possible three momenta Pi 
in unit volume, and prefactor 2 takes into account the initial spin states of the neutron. 
The beta-decay matrix element squared, after summing over spins of final particles and 
averaging over angles, is obtained as G'^cos'^6c{l + 5g\), where G = 1.436 x 10"'^^ erg cm? 
is the weak-coupling constant, 9c is the Cabibbo angle, and qa = —1.261 being the axial 
vector coupling constant. Final electron and proton states must be vacant if the reaction 
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is to occur, and this accounts for the blocking factors I — n2 and 1 — na. The electron- 
capture process (8.4) gives the same energy loss rate as process (8.3), but in neutrinos, 
and therefore the total luminosity per unit volume of the Urea process is twice of Eq. 
(8.8). The integrals may be calculated straightforwardly, since the neutrons, protons, and 
electrons are very degenerate. One thus finds (Boltzmann's constant /c^ = 1) 



Eurca = ^lo^g rrinmpHeiTyef (5.9) 



Here is the threshold factor 0(pe + Pp — Pn), which is +1 if the argument exceeds 0, 
and is otherwise. 

Particle interactions change this result in a number of ways. First, the neutron and 
proton densities of states are determined by effective mass rather than bare masses. Sec- 
ond, the effective weak-interaction matrix elements can be modified by the medium. These 
effects are expected to reduce the luminosity, but probably by less than a factor of 10. 

The temperature dependence of the direct Urea emissivity may easily be understood 
from phase-space considerations. The neutrino or antineutrino momentum is ~ T, and 
thus the phase space available in final states might be a three-dimensional sphere of this 
radius, whose volume is proportional to {T)^- The participating neutrinos, protons, and 
electrons are degenerate. Therefore, for the reaction to occur, they must have energies 
that lie within ~ T of the energies at the Fermi surfaces, and thus each degenerate particle 
contributes a factor ~ T. 

Yet another possibility for the state of matter at high densities is quark matter, in 
which quarks can move around essentially as free particles, rather than being bound 
together as colour singlet entities, such as nucleons and pions. Neutrino emission from 
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such a system was considered by Iwamoto [ |148|| . The basic processes are the quark 
analogue of the nucleon direct Urea processes (8.3) and (8.4), i.e., 



d 



u + e + Ue 



(5.10) 



and 



u + e 



d + z/g. 



(5.11) 



The condition for beta equilibrium is 



(5.12) 



If quarks and electrons are treated as massless noninteracting particles, this condition is 



which is identical to the condition for it to be just possible to conserve momentum for 
excitations near the respective Fermi surfaces. At threshold the momenta of the u quark, 
the d quark, and the electron must be coUinear, but, as Iwamoto pointed out, the weak- 
interaction matrix element for this case vanishes. However, if quark-quark interactions 
are taken into account, the direct Urea process for quarks and electrons that are not 
collinear will be kinematically allowed. To illustrate this effect, consider the case in which 
interactions may be treated perturbatively. To first order in the QCD couphng constant 
a, the quark chemical potentials are given by 



Pf(d) =Pfiu)+pf{e), 



(5.13) 



/ii = [1 + —a]pf{i), i = u, d, 



(5.14) 
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while the electron chemical potential is unchanged. Since a is positive, it is easy to see 
that the conditions for beta equilibrium [Eq. (8.12)] and momentum conservation may be 
satisfied simultaneously. Angles characterizing deviations from coUinearity are typically 
of order a^^^. The calculation of the neutrino and antineutrino emission rates proceeds 
in essentially the same way as for the nucleon process, and the overall results for the 
luminosity is 

Eq = ^^^^^apf{d)pf{u)fx^{Tf. (5.15) 

This has a form similar to the nucleon Urea rate (8.9), but there are some significant 
differences. First, there is a factor a, which refiects the fact mentioned above that the 
weak-interaction matrix element vanishes for collinear relativistic particles, whereas for 
non-relativistic nucleons the corresponding matrix element is essentially independent of 
angle. The second difference is that the quantities Pf{u) and Pf{d) take the place of 
the nucleon masses. Third, the numerical coefficient is different because for quarks the 
angular dependence of the matrix element is important. However, since Pf{u) and Pf{d) 
are expected to be less than m„, and a is less than or of the order of unity, the neutrino 
luminosity from quark matter is expected to be rather less than the characteristic rate for 
the nucleon process. However, it is important to note that for quark matter the electron 
fraction is uncertain. For instance, if u, d, and s quarks may be treated as massless 
and free, the electron fraction vanishes identically. Detailed estimates of the composition 
of quark matter for various models are given by Duncan, Shapiro and Wasserman [p. 49 



and Alcock, Farhi and Olinto |13£]. So far, we have assumed the quark to be massless. 
While this is a good approximation for u and d quarks, but poor for s quarks, which can 
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participate in Urea proeess even in the absence of strong interactions. However, detailed 
calculations show that the energy-loss rate from processes in which s quarks participate is 
less than that from the processes where for u and d quarks are involved (Iwamoto [ |14^ ). 



5.2 Equilibrium neutrino emissivity of quark 
matter 



Iwamoto ||148|| has derived the formula for e using apparently reasonable approxima- 



tions and this formula has been widely used |148, 139, 103, 150| to calculate e for two and 



three flavour quark matter. According to his formula, e is proportional to baryon density 
ub, strong coupling constant and sixth power of temperature T for d quark decay. For 
s quark decay, T dependence of e is same as that for d decay. Furthermore, his results 
imply that electron and quark masses have negligible effect on e and s quark decay (in 
case of three flavour quark matter) which play a rather insigniflcant role. 

In this section, we describe [ [162|| an exact numerical calculation of e and a comparison 



of our results with the Iwamoto formula. Our results show that the Iwamoto formula 
overestimates e by orders of magnitude when Pf{u) +pf{e) —pf{d{s)) is comparable with 
the temperature. For reasonable values of ac and baryon densities, this quantity is much 
larger than the expected temperatures of neutron stars (~ few lOths of MeV) for two 
flavour quark matter, but is comparable with temperature for three flavour quark matter. 
The neutrinos are emitted from the quark matter through reactions 

d ^ u + e~ + Ue 
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u + e' ^ d + Ue 
s u + e~ + Ve 

M + e~ — >• s + z/e- (5.16) 

The equilibrium constitution of the quark matter is determined by its baryon density ub, 
charge neutrahty conditions and weak interactions given in Eq. (8.16). Thus, for two 
flavour quark matter, 

IJ'd = f^U + f^eifJ'Ue = f^Pe = 0) 

2n„ - nrf - 3ne = 
nB = {nu + nd)/3 (5.17) 

and for three flavour quark matter 

IJ'd = IJ'u+ f^eif^ue = f^Pe = 0) 
/id = l^s 

2nu - Ud - Us - Sue = 
Ub = {nu + nd + ns)/3. (5.18) 

The number density of species i is rii = g .p^jr{i) / {6tt'^) with the degeneracy factor gi being 



two for electron and six for quarks. For electrons fie = yP}{^) + ^'^^ quarks we 



use 
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f^,= [^ + ^{l--ln{x + v))]Pf, (5.19) 

where x = Pf{q)/mg and r] = Vl + x^, f^q being the quark mass. For massless quarks, 
Eq. (8.19) reduces to 
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So- 

^Ji, = {l + ^)Pf{q). (5.20) 

The neutrino emissivity e for reactions involving d{s) quarks is calculated by using the 
reactions in Eq. (8.16). In terms of the reaction rates of these equations, we get ( 

n = c = I), 

{Pd{s)-Pu){Pu-Pe) 

EuEd{s)Ee 



^d{s) = Ad(s) / d pd{s)d Pud Ped p, 



xS'^iPdis) - Pu - Pe - Pu)n{pd{s))[^ - n{pu)][l - n{pe)], (5.21) 

where Pi = {Ei, pi) are the four momenta of the particles, n(pj) = are the 

Fermi distribution functions and 

24G2cos2^, 

^^-^^ (5-22) 

For degenerate particles, {[3pf{i) ^1), Iwamoto has evaluated the integrals in Eq. (8.21) 
using certain reasonable approximations and obtained the simple expressions for and 



as given below 148 



914 

Q = —G"^ cos'^ ecacPf{d)pf{u)pf{e)T^ 

C sin^ e,a,fisPf{u)pf{e)T^ (5.24) 
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The approximations involved in obtaining these formulas are 



1. neglect of neutrino momentum in momentum conservation, 

2. replacing the matrix elements by some angle averaged value, 
and 
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3. decoupling momentum and angle integrals. 

The expressions for neutrino emissivity as obtained by Iwamoto have been used widely. 
The temperature dependence of emissivity as obtained by Iwamoto has a physical expla- 
nation. Each degenerate fermion gives one power of T from the phase space integral ( 
(Ppi — > pf{if'dEi (X T ). Thus one gets from quarks and electrons. Phase space 
integral for the neutrino gives d^pi, oc {E^)dEi, oc T^. Energy conserving 5— function 
gives one which is cancelled by one ex T factor coming from matrix ele- 
ment. So finally, one gets e oc T^. This argument, however, ignores the fact that 
Apd {Aps) = pf{u) + pf{e) — pf{d) (p/(s)), which is related to the angle between 
Pd, Pu and Pe could be small and comparable to T. We shall demonstrate below that 
precisely in this region that the Iwamoto formula fails. 

Before discussing the causes of the shortcoming of Iwamoto formula, let us first com- 
pare our results with the Iwamoto formula and try to find out the specific cases where the 
deviation is more pronounced. In Figs. 8. 1-8.3, we have plotted e vs T for two-fiavour d 
decay, three-flavour d decay and s decay respectively. For two-flavour d decay our results 
(e^) are in good agreement with the emissivity calculated using Iwamoto formula (e^/). In 
Fig. 8.1 curves (a) and (b) are eai and respectively, for etc = 0.1 and ub = 0.4. (c) and 
(d) corresponds to the same but for ac— 0.1 and ub— 1.4. It is evident from the flgure 
that agreement of Iwamoto results with our calculation is better for higher densities and 
lower temperatures. Also is consistently smaller than eai, the Iwamoto result, in the 
range of temperatures considered. Corresponding Fermi momenta of quarks and electrons 
are given in Table 8.1. It is to be noted that all the momenta are much larger than the 
temperature. 
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Fig. 8.1: Two flavour rf-decay for = 0.1; (a) Iwamoto results for Hb = 0.4/m ^, (b) 
Our results for ub = {^Pd = 6.78), (c) Iwamoto results for hb = lAfm^^, (d) 

Our results for ub = lAfm'^ {Apd = 10.29). 
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Fig. 8.2: Three flavour rf-decay for etc = 0.1 and s quark mass is 150 MeV; (a) Our 
results for ub = 1.4/m~^, (b) Iwamoto results for ub = 1.4/m~^ {^Pd = 0.067), (c) Our 
results for = 0.4/m~^, (d) Iwamoto results for = 0.4/m~^ {^Pd — 0.39). 
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Fig. 8.3: Three flavour s-decay for = 0.1 and s quark mass is 150 MeV; (a) Our 
results for ub = 1.4/m~^, (b) Iwamoto results for ub = 1.4/m^^ {^Ps = 1-613), (c) Our 
results for ub — 0.4/m~^, (d) Iwamoto results for ub — OAfm'^ {^Pd — 9.719). 
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TABLE 8.1 



BARYON NUMBER DENSITY n^, FERMI MOMENTA OF m-QUARK pf{u), 
d-QUARK pf{d), AND ELECTRON pf{e) FOR DIFFERENT ac, 
WHERE Apd =pf{u)+pf{e)-pf{d). 







Pf{u) 


Pf{d) 


P/(e) 








(MeV) 


(MeV) 


(MeV) 


(MeV) 




0.60 


357.80 


449.20 


99.15 


7.75 


0.1 


1.00 


424.22 


532.52 


117.56 


9.20 




1.40 


474.57 


595.79 


131.51 


10.29 




0.60 


357.71 


449.25 


95.43 


3.89 


.05 


1.00 


424.11 


532.65 


113.15 


4.61 




1.40 


474.45 


595.87 


126.57 


5.15 



Fig. 8.2 shows the e^, for three-flavour quark matter. It shows that e^i is 2 -3 orders 
of magnitude higher compared to our results. Here contrary to the two flavour case, the 
difference becomes more pronounced at higher densities. Fig. 8.3 shows the variation of 
Es with temperature. Here again, it is clear that is quite different from e^/ but this 
difference is less compared to that between and e^/. For all the cases, the difference 
between our results and those using Iwamoto formula increases at higher temperatures. 
The Fermi momenta of quarks and electron for three- flavour case are given in Table 
8.2. The study of all the figures and tables above reveals that the cases where Iwamoto 
formula agrees reasonably well with our results, Ap^ (or Aps ) is much larger than the 
temperature. On the other hand, when this difference is smaller or comparable with the 
temperature, the Iwamoto formula overestimates the exact result by order of magnitude. 
In addition to these. Fig. 8.3 also shows that our results are about a factor of 2.5 lower 
than the Iwamoto results even at lower temperatures. We have found that this difference 
comes from the approximation involved in the calculation of matrix element. 
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TABLE 8.2 

BARYON NUMBER DENSITY n^, FERMI MOMENTA OF ?i-QUARK pf{u), d-QUARK 
Pf{d), s-QUARK pf{s) AND ELECTRON pf{e) FOR DIFFERENT m, AND 
DIFFERENT a^, WHERE Apd = pf{u)+pf{e) -pf{d) 
AND Aps=pfiu) +Pf{e) -pf{d) 



nis 


ac 




Pf{u) 


Pf{d) 




P/(e) 




Aps 


{MeV) 






(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


{MeV) 






0.60 


356.99 


360.06 


353.86 


3.33 


0.26 


6.46 




0.1 


1.00 


423.26 


424.81 


421.69 


1.69 


0.14 


3.26 


150.0 




1.40 


473.49 


474.27 


472.72 


0.84 


0.06 


1.61 






0.60 


356.99 


365.89 


347.62 


9.28 


0.38 


18.65 




0.05 


1.00 


423.26 


430.29 


415.98 


7.33 


0.30 


14.61 






1.40 


473.49 


479.49 


467.34 


6.25 


0.25 


12.40 






0.60 


356.99 


365.93 


347.58 


9.70 


0.76 


19.11 




0.1 


1.00 


423.26 


429.10 


417.25 


6.34 


0.50 


12.35 


200.0 




1.40 


473.49 


477.66 


469.25 


4.52 


0.35 


8.76 






0.60 


356.99 


374.26 


337.85 


18.00 


0.73 


37.14 




0.05 


1.00 


423.26 


437.14 


408.39 


14.47 


0.59 


29.34 






1.40 


473.49 


485.45 


460.90 


12.46 


0.50 


25.05 
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Furthermore, Table 8.2. shows that for three flavour case electron chemical potential 
(which is same as p/(e) for massless electrons ) becomes small (< 1 MeV ) for some values 
of ac, ub and m^. In these cases, electrons are no longer degenerate. Clearly, for such 
cases the Iwamoto formula is not applicable. This point is missed in earlier calculations. 

Our results are expected to have important implications for neutrino emissivity and 
quark star cooling rates because all the earlier calculations have used the Iwamoto formula 
and predicted large quark star cooling rates in comparison with the neutron star cooling 
rates for temperatures less than 1 MeV. Our results show that particularly, for three- 
flavour quark matter, the calculated emissivity is at least two orders of magnitude smaller 
than the one given by Iwamoto formula and therefore, the three-flavour quark star cooling 
rates are that much smaller. Hence, it is necessary to understand why Iwamoto formula 
fails. 

To investigate the shortcoming of the Iwamoto formula, we consider the integral 



Here, we have replaced the neutrino emission rate by unity and therefore / is essentially 
the phase space integral. Following the reasoning of Iwamoto, this integral should be 
proportional to T^. Choosing the coordinate axes such that pd is along z-axis and Pu is in 
X — z plane and using the three-momentum 5— function to perform electron and li-quark 
angle integrations, we get 




x^^{Pd{s) - Pu - Pe - Pu)n{pd{s))[l - n{pu)][l - n{pe)]. 



(5.25) 
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0^ 



xS{ed - eu - ee - e„)n{pd)[l - n{pu)][l - n{pe)], 



(5.26) 



where Xi, — cos 61, and Xu — cos 6u is determined by solving 



Pu^ 



u 



= Pd-PuXu - 



\pI-pI{i-xI)-pI{i-xI) 



2p^p,.J{l-xl){l-xl)cos4>,Y'\ 



(5.27) 



The integral in Eq. (8.26) above is restricted to the momenta \pi — Pf{i)\ few times T 
due to Fermi distribution functions and the energy (5— function. Now, if we neglect the 
neutrino momentum in the 5— functions, we get, Xu = iPd'''' Pu~ pI) f^PdPu and the factor 
in the square brackets of Eq. (8.26) becomes 1/pdPuPe- 
Two points should be noted at this stage. 

1. Generally, Xu is close to unity, so that 1 — is small. But, if Apd is of the order of T, 



1 — xlpd can be comparable with T and p^ and therefore p^ cannot be neglected 



in the momentum 5— functions. Particularly, the denominator in the square bracket 



of Eq. (8.26) cannot be approximated bypePuPdyl ~ ^u- Thus, if p^y 1 ~ < Pu, 
one would get a power of T from the denominator and / will not be proportional to 

2. Secondly, the momenta may differ from the corresponding Fermi momenta by few 
times T in the integral. When Apd ~ T, there are regions in pdPuPe—spa.ce where 
Xu > 1 and the rest of the integrand is not small. Clearly, these regions must be 
excluded from the integration as these values of x^ are unphysical. If one does not 
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put this restriction, as is done when one factorises angle and momentum integrals, 
the phase space integral will be overestimated (and wrong) when Ap^^ ~ T. 



The above discussion clearly shows why the integral in Eq. (8.26) should not be pro- 
portional to when Ap^ ~ T. In order to demonstrate this point, we have calculated 
the integral in Eq. (8.26) numerically and compared with the approximation where the 
neutrino momenta are neglected and the restriction imposed by Xu condition is not im- 
posed. The calculation is done for = 0.1 and for two- flavour case. The results are 
shown in Fig. 8.4. In this figure, we also show the result for a case where the electron 
mass is taken to be 25 MeV. This is of course unphysical, but by adjusting the elec- 
tron mass we can reduce Ap^. The figure clearly shows that the approximate value of 
I is proportional to where as the exact integral is smaller than the approximate one 
at large T. Further more, for 25 MeV electron mass, the departure from sets in at 
smaller value of the temperature. This clearly shows that the departure is dependent on 
the value of Ap^. Here we would like to mention that for some values of etc and , P/(e) 
is small and is of the order of T. This implies that electrons are no longer degenerate and 
deviation from the Iwamoto result is most pronounced. 

In Eq. (8.26), we have dropped the matrix element of the weak interaction in the 
emissivity calculation (Eq. (8.21)). So, the discussion of preceeding paragraphs apply to 
the emissivity calculation as well. Therefore, it is now clear why the Iwamoto formula 
fails, when Apd (or Aps in case of weak interactions involving strange quarks) is close to 
the temperature of the quark matter. 



Similar approximations have been used by other authors [|164| , |165|] . In Ref. ||164 
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Fig. 8.4: Two flavour phase space integrals for = 0.1 (a) Without restriction on cos^„ 
for both electron mass me=0.0 and 25 MeV , (b) Exact integral for me=0.0, (c) Exact 
integral for me—25 MeV. 
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Fig. 8.5: is plotted against x where x — . ^ . The fitted function 

1 + ax + bx'^ + cx^ where a = -2.5, 6 = 100. and c = 30. 



90 



Burrows has calculated the neutrino emissivity for non-interacting quark matter from re- 
verse beta decay. The exponent of T is 7 instead of 6 (Iwamoto), because of the partial 
restriction of the electron's phase space. Duncan et al. | |165| | computed the emissivity 



from reverse and direct beta decay for both interacting and non-interacting quark matter. 
They have reproduced the results of Burrows (T'') [|164|] for non-interacting reverse beta 



decay and Iwamoto (T ) ||148|| for interacting quark matter. In both these works, the 
effect of the finite neutrino momentum is included but rest of the calculation follows the 
approximation scheme of Iwamoto. Hence, their temperature dependence of emissivity 
formula is different from ours. Here, we would like to note that the departure from 
dependence of the emissivity essentially arises from the careful phase space integration. 
Since, similar approximation scheme is used to obtain the neutrino emissivity of neu- 
tron matter, it is possible that the emissivity calculated for neutron matter may also be 
overestimated when x is large. 

Since the departure from the Iwamoto formula arises from the fact that T/Apd (or 
T/ Aps ) is not small, it may be possible to fit the numerically calculated e with a function 
of the form ej/f{x), where x = T / Ap^ {T/Aps for strange sector). The function f{x) 
should be such that for small values of x it should approach unity. Choosing /(x) = 
1 + ax + + cx^, we have fitted the calculated e for a number of values oi ub-, cxc 
and nis for both s and d decay and obtained the values of a, h and c. The quality of fit is 
shown in Fig. 8.5 (since, for s- decay, as mentioned above, there is difference of a factor of 
2.5 in Iwamoto and our results even at lower temperatures, the data points for s- decay, 
in Fig. 8.5., have been scaled accordingly). The values of a, h and c are —2.5, 100. and 
30. respectively. It is clear from Fig. 8.5 that for x ^ 0, Iwamoto results approache to 
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our values. Hence, our fitting is valid for any values of x. 
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5.3 Non-equilibrium neutrino emissivity of 
quark matter 



The neutrino emission is a dominating mechanism of coohng of quark or neutron 
stars with the internal temperature exceeding ~ 10^ K. The quark stars cool much faster 



than the neutron stars. The neutrino emissivity formula ||148| , |162|| of the quark stars are 
obtained under the assumption that the system is in equilibrium with respect to weak 
interactions. However, as long as matter is strongly degenerate, the number density of 
neutrinos is small and their presence in the matter can be neglected. Thus, the beta 
equilibrium (chemical equilibrium) stands as a good approximation, which implies the 
equality of anti- neutrino and neutrino emissivity. Of course, one can work even though 
the chemical equilibrium is not satisfied. Radial pulsations of quark stars have periods 
~ milliseconds , much shorter than the time scale of the beta reaction at T < lO^^K ~ 
1 MeV. Local compression as well as the rarefaction are very much important for non- 
equilibrium beta reactions and depend on the difference between the chemical potentials 
of u, d{s) quarks and e electron, i.e., Sfi = Hd{s) — l^u — l^e- The non-equilibrium beta 
reactions, induced by the radial pulsations of neutron and quark stars, were studied by 
a number of authors ||130| , |153|| . Recently, Madsen ||155|| and Sawyer ||151|| calculated the 



bulk viscosity and damping rates by invoking the non-equilibrium condition in the quark 



star matter. However, later, Haensel ||166|| calculated the neutrino emissivity of non- 



equilibrium neutron star matter, where he found that the emissivity depends on the non- 
equilibrium conditions such as the difference in the chemical potential of neutron, proton 
and electron = yUn — — /^e)- Moreover, he showed that with increase of 5^ {—5^') 
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the emissivity increases (decreases). So, it is interesting to see the non-equihbrium effects 
in the quark star matter. 

In this section, we consider the rapid compression of the hquid interior in quark stars, 
which has two/three component quark matter. Such a situation is expected to occur 



during the gravitational collapse of the neutron stars or the quark stars in black hole ||167 
which would take place when the mass of an accreting neutron star exceeds the maximum 
allowable mass for the equilibrium configurations. Thus, a significant deviation from the 
chemical equilibrium is to be expected because of the shrinking of the stellar radius, 
which implies the monotonic increase of the average density in collapsing stars. Here 



we concentrate ||168|| on the characteristic features of the beta non-equilibrium neutrino 
spectra and compared with that of the beta equilibrium neutrino spectra. 

Let us consider the neutrino emissivity of quark matter having two/three flavour de- 
grees of freedom. Each of the constituents separately is a Fermi liquid in thermodynamic 
equilibrium. 

For the degenerate two flavour quark matter, the simplest neutrino processes are the 
direct beta decay reactions 

direct [3- : d ^ u + e~ + Ve 

direct (3+ : u + e~ d + v^- (5.28) 

The distinction between /3_ and /3+ processes will be discussed later on. 
Moreover, the charge neutrality process of the two flavour quark matter is 

2nu — rid — 3Tie = (5.29) 
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and the baryon density is defined as = (n^ + no)/ 3 . Since the matter is strongly 
degenerate, we write the Fermi momentum ppi — (67r^nj/g')^/^ by neglecting the thermal 
corrections, being the number density of i-th particle and the degeneracy factor g is 
6 for quark and 2 for electron. However, the non-equilibrium condition for the above 
reaction is 

/^d - - /^e = (^/^(/^l^e = /^Pe = 0). (5.30) 

Here, we approximate the positive 5/1 for the direct /3_ reaction (anti-neutrino emission) 
and that of the negative S/j, for the direct /3+ reaction (neutrino emission). 
Similarly, for the three flavour degenerate quark matter, the simplest neutrino process are 
the direct beta decay reaction, which occurs through d and s quarks. In addition to Eq. 
(8.28), one has also the following reactions : 

direct /3_ : s — > m + e~ + i/g 

direct P+ : u + e~ ^ s + 1/^ (5.31) 
The equations for the charge neutrality and the baryon density are respectively 

2n„ — rid — ris — 3^e = (5.32) 

and 

ns = (riu + nd + ns)/3, (5.33) 

and due to the non-equilibrium condition, one can write 

l^d- l^s = (5.34) 
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Moreover, here we again approximate the positive for the direct (3- reaction (anti- 
neutrino emission) and the negative 5/i for that of the direct /5+ reaction (neutrino emis- 



sion . 



Thus, in these approximations for both the two and the three flavour quark matter, the 
neutrino emissivity is almost same, i.e., ep^{~6fi) ^ 6u^{6fi). In the chemical equilibrium, 
= 0, which implies Sp^ ~ e^^. 

The energy momentum relation of s quark is approximated by Eq. (8.19) and that of 
u and d quarks by Eq. (8.20), i.e., /i is replaced by E and pp is replaced by p, in above 
Eqs. (8.19-8.20). This is reasonable since only the energies and the momenta close to the 
Fermi surface will contribute to the matrix element. The neutrino emissivity 6d(s) for d{s) 



decay | 148 | is given by Eq. (8.21). The momentum delta function in Eq. (8.21) is used 
to integrate over the neutrino momentum p^. For the massless electrons, = Pe and the 
energy delta function is used to perform the integral over p^. This gives 

^ _ (-gd(s) - EuY - {Pdjs) -Puf - 2pd(^s)Pu{l - COsOu) 
{2puCOs6ue + 2pd{s)COs6e + 2Ed(s) - 2Eu) 

where 9u, Oe and Que are the angles between d and u, d and e and u and e respectively. 

The chemical potentials and Fermi momenta thus obtained by employing the above 
non-equilibrium and charge neutrality conditions for two and three flavour matter are 
substituted in neutrino emissivity expression [Eq. 8.21] which subsequently yield 

3^ . . ^3^ ^2^n iVd{s)-Vv)ij>u-Ve) 
EuEd{s)Ee 

xn(pd(.))[l - n{pu)][l - n{p^)\. (5.36) 



ed(s) = Ad{s) / d pd{s)d PuPed^ 



In the present case, the above emissivity has been evaluated numerically. The integral 
is five dimensional as all the angles are measured with respect to d{s) quark. Throughout 

96 



the calculation, we have taken u, d, e, and u masses to be zero, s quark mass and 
strong coupling constant to be 150 MeV and 0.1 respectively. 

Recently [|162|| , the neutrino emissivity of the degenerate quark matter was calculated 
exactly with the chemical equilibrium {6fi = 0.0), which has been discussed in the previous 
section. It has been found that the neutrino emissivity results are in qualitative agreement 



with that of Iwamoto ||148|| for the two flavour quark matter, whereas for the three flavour 
quark matter Iwamoto's result overestimates the numerical values by nearly 2 orders of 
magnitude or more for d decay and agrees with that of the s-decay results within a factor 
of 3-4. Also, it was pointed out that the dependence of the temperature and the density 
on the neutrino emissivity is quite different from Iwamoto's result and is sensitive to the 
s quark mass. 

Here, we use the same method [ [L62|| to calculate the neutrino emissivity by employing 
the non-equilibrium condition. For the two values of temperatures as well as baryon 
densities and different values of 5/1, the numerical values of the neutrino emissivity for 
two and three flavour quark matter are quoted in Table 8.3 and 8.4 respectively. The 
variation of the corresponding Fermi momenta with respect to the increase and decrease 
of are presented in Table 8.5 and 8.6. We found here that the dependence of emissivity 
on ub and T is rather different from that of the equilibrium condition {6fi = 0). 

For the two flavour quark matter, the anti-neutrino emissivity increases while that of 
the neutrino emissivity decreases on increasing 6fi (same as increasing —6fi in the case 
of neutrino). For a constant temperature of T = 0.4 MeV, the anti-neutrino emissivity 
depends on its exponent 

2^3.02 2^1.47 that of the baryon density as 

Ep OC n^^ to n}^"^ for 5/i varies from 1. to 2. MeV (0.005 to O.Ol/m^i) In a similar 
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TABLE 8.3 

EMISSIVITY FOR TWO FLAVOUR QUARK MATTER. HERE e^(5/x) AND 
e'^{-5n) ARE THE EMISSIVITIES FOR d DECAY. 



{MeV) {fm"^) {erg / cm? / s) {erg / cm? / s) fm~^ 

(L8 7.64x10+^^ 1.20x10+^^ 0005" 

1.0 9.59x10+^^ 1.54xl0+-'^^ 

M 7.22x10+^'^ 4.50x10+^^ 0.0075 

1.0 9.02x10+^6 5.80x10+^^ 

0.4 (TS 3.75x10+^^ 1.36x10+^^ (ToF 

1.0 4.52x10+^^ 1.79x10+^^ 

M 4.50x10+^^ 4.50x10+^^ OW 

1.0 5.70x10+^^ 5.70x10+^^ 

(18 6^67x10+28 2.33x10+^^ 0.005 

1.0 8.41x10+28 3.02x10+-^^ 

(l8 1.12xl0+^« TJOxW^ 0.0075 

1.0 1.42x10+28 9.29x10+-^^ 

0.5 (l8 1.58x10+2'^ 1.86x10+^2 q^qF 

1.0 1.97x10+27 2.49xl0+-'^2 

M 1.61x10+^^ 1.61x10+^*^ 000" 

1.0 2.05x10+^0 2.05x10+=^° 
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TABLE 8.4 

EMISSIVITY FOR THREE FLAVOUR QUARK MATTER. HERE e^(5/x) AND e^(-(5/i) 
ARE THE EMISSIVITIES FOR d DECAY AND e^(5^) AND ei(-5//) 
ARE THE EMISSIVITIES FOR s DECAY. 



T ub 




ei{-5n) el{5^. 




(MeV) ifm-^) 


(erg/cm^/s) 


[er g 1 cm? j s) {erg /cm 


I'^/s) {erg/cw? / s) MeV 



(TS L54xT0+23 L46xT0+23 6.47x10+^^ 8.52x10+^^ LOO 
1.0 8.31x10+23 1.78x10+23 4.20x10+26 6.51x10+23 



(TS 2.55x10+^^ 2.01x10+^3 ^MxW^ 1.26x10+^^ 

1.0 1.51x10+2^ 6.51x10+22 3.00x10+27 9.87x10+2^ 

0.2 08 3.75x10+^* 7.76x10+^^ 2.43xlO+2« L23><T0+^" 3l)0 

1.0 2.53x10+2^ 1.82x10+22 1.07x10+2^ 9.75x10+^^ 

08 8.72x10+23 8.72x10+23 3.43x10+2^ 3.43x10+2^ OOO 

1.0 l.lOxlO+2-^ l.lOxlO+2-^ 2.19x10+2^'' 2.19x10+2^^' 

08 3.02x10+25 1.20x10+25 6.18x10+2^ 2.64xl0+2« LOO 

1.0 1.94x10+25 6.93x10+2^ 3.20x10+2^ 1.58x10+2^ 

08 450x10+25 7.08x10+24 2.07x10+2* 4.08x10+25 2m 

1.0 3.06x10+2^ 3.75x10+2^ 9.73x10+2^ 2.49x10+2^ 

0.4 08 6.49x10+25 3.93x10+^^ 5.49xl0+2« 5.45x10+2^ SlX) 

1.0 4.60x10+25 1.80x10+2^ 2.33x10+2^ 3.41x10+2^ 

08 1.93x10+^5 1.93x10+25 1.43x10+2^ 1.43x10+2^ OOO 

1.0 1.19x10+25 1.19x10+25 8.15x10+26 8.15x10+26 
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way, the neutrino emissivity is proportional to the temperature as well as the baryon 
density exponents (for the constant temperature of T = 0.4 MeV as in the earlier case) 
as oc T9-^2 ^.^^ ^im ^0.84 ^ -2. MeV (0.005 to 0.01 fm-1). 

But for the equilibrium case {6fi = 0), ep^ oc T^-^-^-^ and oc n^^^^ as has been observed 
earlier [|162|. 



In addition, for the three flavour quark matter, the behaviour of the anti-neutrino as 
well as the neutrino emissivity is qualitatively same due to the variation in 6fi{-6fi in the 
case of neutrino), as is seen above in the case of two flavour quark matter. But unlike 
the two flavour quark matter, the anti-neutrino and the neutrino emissivity is inversely 
proportional to the baryon density. The variation in the exponents here (for T = 0.4 MeV 
in both cases) is, oc T^-^^ to T^-^^ and n]^^'^^ to n^'^^ for the d decay, and oc T^'^^ 
to T^'^^; and n]^''^^ to n]^'^^ for the s decay when one changes 611 from 1. to 2. MeV . 
Similarly, oc T^-^s to T^-^^ and n^j^-^s to nf -^^ for the d decay, and oc T^-^^ to T^^-^^; 
and n'^ '^^ to n^'^"^ for the s decay due to variation in 5/1= -1. to -2. MeV . But the 
results in the equilibrium case (J/i = 0) which has already been shown in the literature 

1 1621 is oc T^-86-4.47. ^ ^^2.16 ^ ^^^^^ ^ rj,^m-^.Zfi ^ ^-2.53 ^ 



decay with = 150 MeV and Oc = 0.1. 

Therefore, it is obvious from the above numerical results that the exponent of the tem- 
perature decreases but that of the baryon density increases in the anti-neutrino emissivity 
case compared to the equilibrium condition. But for the case of neutrino emissivity, the 
exponent of the temperature increases and the baryon density exponent decreases. Thus 
it is clear that the changes on the neutrino emissivity as well as the anti-neutrino emis- 
sivity is due to the non-equilibrium condition where 7^ 0. The change in the emissivity 
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{£p{e,y)) due to the variation in 5ii{—5ii) for T = 0.4 MeV is shown in Table 8.3 and 8.4 
for two flavour and three flavour case, which has also been extended to T = 0.5 MeV 
(Table 8.3 and 8.4). Moreover, it is seen from Table 8.5 and 8.6 that the difference in 
the chemical potential affects the density distribution of particles for both two and three 
flavour quark matter case. So, when increases from 0.0 to 0.01 fm~^, the Fermi 

momentum of u quark and e increase (decrease) , and that of d quark decreases (increases) 
with respect to the chemical equilibrium Fermi momenta oi u, d and e for two flavour 
quark matter. But for the three flavour quark matter case (for d and s decay), the chem- 
ical equilibrium Fermi momentum of u quark remains the same, and that of d quark and 
e decrease (increase) and s increases (decreases) as one changes 5iJt,{—5ij) from 0.0 to 2.0 
MeV . Hence, in comparison to the equilibrium case, the neutrino emissivity decreases 
monotonically, whereas that of the anti-neutrino emissivity increases monotonically for 
both the two and the three flavour quark matter. 
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TABLE 8.5 

BARYON NUMBER DENSITY hb, FERMI MOMENTA OF li-QUARK 
Pf{u), d-QUARK pF{d) AND ELECTRON pj.(e). 



5/1 


Ub 


Pf{u) 


PF{d) 


PF{e) 






(MeV) 


(MeV) 


(MeV) 




0.8 


393.81 


494.40 


109.13 


0.0 


1.0 


424.22 


532.58 


117.56 




0.8 


393.89 


494.35 


112.16 


0.0075 


1.0 


424.30 


532.53 


120.59 




0.8 


393.92 


494.34 


113.17 


0.01 


1.0 


424.33 


532.51 


121.60 




0.8 


393.74 


494.45 


106.09 


-.0075 


1.0 


424.14 


532.63 


114.52 




0.8 


393.71 


494.47 


105.08 


-.01 


1.0 


424.12 


532.64 


113.50 



TABLE 8.6 

BARYON NUMBER DENSITY n^, FERMI MOMENTA OF m-QUARK pf{u), d-QUARK 
PF{d), 5-QUARK pf{s) AND ELECTRON p^(e) FOR m, = 150 MeV AND = 0.1. 





Ub 


Pf{u) 


PF{d) 




PFie) 


(MeV) 




{MeV) 


{MeV) 


(MeV) 


(MeV) 




0.8 


392.92 


395.08 


390.73 


2.35 


0.0 


0.1 


123.26 


121.81 


121.69 


1.69 




0.8 


392.92 


394.62 


391.20 


1.85 


1.0 


0.1 


423.26 


424.35 


422.16 


1.18 




0.8 


392.92 


394.15 


391.68 


1.34 


2.0 


0.1 


423.26 


423.88 


422.63 


0.68 




0.8 


392.92 


395.55 


390.25 


2.86 


-1.0 


0.1 


423.26 


425.28 


421.21 


2.19 




0.8 


392.92 


396.01 


389.77 


3.36 


-2.0 


0.1 


423.26 


425.74 


420.74 


2.70 
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Chapter 6 



Summary and Conclusion 



In this thesis, we have investigated the behaviour of matter at very high density 
using a relativistic Lagrangian description. The Lagrangian chosen by us corresponds to 
the chiral sigma model. This approach is considered to be a "good" low energy hmit 
of quantum chromo dynamics. Although there have been a few previous calculations 
along this line, a detailed and consistent field theoretical approach has been lacking. 
The calculations presented in this thesis are aimed at such a detailed study. We have 
extended these calculations for the case of finite temperatures (< 15 MeV). The results 
are expected to find application in stellar collapse calculations. In addition, we have 
dealt with the following subjects : (1) Phase transition to quark matter and the possible 
formation of strangelets at high densities and (2) astrophysical applications of our results, 
to (a) structure and radial oscillation of nonrotating neutron stars and (b) the neutrino 
emissivity of quark matter with an improved calculation of phase space integrals involved. 
The highlights and main results of this thesis can be summarized as follows : 

1. The energy per nucleon of cold nuclear matter {ksT — 0), derived by us using 
chiral sigma model, is in good agreement with the preliminary estimates inferred 
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from heavy-ion collision data |0] in the density range between one to four times 



the nuclear saturation density (ris). 

2. For a system of high density nuclear matter, based on the chiral sigma model, we find 
that a strict first order phase transition to [u, d, s) quark matter is not favoured. 
This does not, of course, preclude a phase transition of second order. However, we 
have not investigated the latter problem. 

3. The mass formulae for finite lumps of strange quark matter with u, d and s quarks 
and non-strange quark matter {u and d) are derived in a non-relativistic approach, 
taking into account the finite size effects such as surface and curvature. We find 
that there is a good possibility for the formation of metastable strangelets of large 
mass detectable in experiment. This is important since the detection of strangelets 
may be the most unambiguous way to confirm the formation of quark-gluon plasma 
in heavy ion collision experiment. 

4. The maximum mass for stable neutron stars predicted by our equation of state 
for (n, p, e) matter is 2.59 times the solar mass. The corresponding radius (i?), 
crustal length (A) and surface red shift ratio (a) are 14.03 fcm, 1.0 km and 0.674 
respectively. The maximum moment of inertia is 4.79 x lO'^^ g cra^ . These suggest 
that our equation of state for neutron star matter is comparatively "stiff' . This is 
reflected in the value of the maximum mass of neutron stars, which is the largest for 
the present model as compared to other available field theoretical equation of state 
models. It may be mentioned here that observational evidence in favour of a stiff 
equation of state comes from the identification by Triimper et al. |Q of the 35 day 
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cycle of the pulsating X-ray source Her X-1 as originating in free precession of the 



rotating neutron star (Pines for a discussion). For a 1.4 times the solar mass 
neutron star configuration, we get : R = 14.77 km, I = 2.15 x 10'^^ g cm?, A=3.0 
km and the red shift ratio (at the surface) a = 0.85. The corresponding central 
density is 4.06 x 10^^ g cm~^. 

5. The neutrino emissivity from two and three flavour quark matter is numerically 



calculated and compared with the result given by Iwamoto [p.48|] . We find that the 
emissivity is smaller than Iwamoto's result by about two orders of magnitude when 
Pf{u) +pf{e) —pf{d{s)) is comparable to the temperature. We attribute this to the 
severe restriction imposed by momentum conservation on the phase space integral. 
An alternative formula for the neutrino emissivity, which is valid when the quarks 
and electrons are degenerate and any values of Pf{u) + pf{e) — pf{d{s)) is obtained 
by us. 
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